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(57) Abstract 

A method of determining or predicting a value Pi of a property (e.g., octane number) of a material X or a property of a product of 
a process from said material or yield of said process for example a blending, separation, or chemical (e.g., polymerisation) process, which 
method comprises measuring the absorption Du of said material at more than one wavelength in the region 600-2600 nm, comparing the 
said signals indicative f said absorptions or a derivative thereof with signals indicative of absorptions Dim or derivatives thereof at the 
same wavelength for a number of standards S in a bank for which the said property or yield P is known, and choosing from the bank at 
least one standard Sm with property Pn said standard having the smallest average value of the absolute difference at each wavelength i 
between signal for the material and the signal for the standard Sm to obtain Fx. with averaging of said properties or yields Pm when 
more than one standard Sm is chosen. 
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PROPERTY DETERMINATION 

This invention relates to a method of determining or predicting by near infra 
red (NIR) spectroscopy properties of feeds or products and/or yields in physical or 
chemical processes or separations, in particular involving hydrocarbons, especially in 
hydrocarbon refineries or for lubricant uses, or chemical processes including 
5 polymerisation. 

NIR spectroscopy has many advantages over other methods of analysis e.g. in 
refineries and can cover a large number of repetitive applications accurately, quickly 
and on line The NIR region between 800 and 2500nm contains the totality of 
molecular information in the form of combinations and overtones from polyatomic 

1 0 vibrations, but Mathematical techniques are needed to exploit this information and to 
calculate the desired parameters EP-A-28525 i . 304232, 305090. the disclosure of 
which is hereby incorporated by reference, describe the use of NIR for determining 
octane number of a product, or determining yields and/or properties of a product of a 
chemical process in a refinery or separation process from analysis on the feeds to that 

1 5 process, and yields and/or propenies of a product of a blending operation again from 
analysis on the feed thereto. 

At present, numerical methods described for modelling physicochemical 
properties based on NIR spectra are all of a correlative nature and involve relations of 
a regressional character between the property(ies) studied Among these multivariable 

20 analyses are multilinear regression (MLR), Principle Component Regression (PLR). 

Canonic regression, and regression by Partial Least Squares (PLS). In all cases there is 
sought between the property and the NIR spectrum a relation which may be linear but 
is usually quadratic or of higher algebraic form involving regression coefficients 
^applied to each absorption. The establishment of any regression requires a progressive 

25 calibration, as the approach is empirical and not supported by a theory. 
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These techniques have disadvantages, the chief of which is the need for 
establishing a strong correlation between the spectainn and the property, and their 
difficulty in dealing with positive or negative synergy between components 
contributing to that property For example for determining chemical composition e.g. 
5 LINA (linear, isoparaffm, Naphthenic. Aromatics) in a hydrocarbon feed to a catalyst 
reformer, a PLS technique based on the NIR spectra has been described for use. The 
model works well on the calibration set but the response of the models when pure 
hydrocarbons are added e g cyclohexane is not satisfactory, as the model predicts 
changes in isoparaffms and naphthenes the reverse of that found experimentally 
1 0 Funhermore there are other practical difficulties, mainly in the need to identify samples 
of families having the same kind of relation between the spectra and the properties to 
be modelled Thus the model may be limited especially with a non linear relation 
between spectrum and property Especially when at the edges of the available data the 
accuracy of the model diminishes The stability of the model is also a problem, as is 
1 5 the need when adding new standards to do laborious revisions to give the new model, 
especially when adjusting to a new feedstock for a process, thus testing 6 properties on 
4 products leaving a distillation unit requires 24 models, each of which has to be 
changed for each change of the feed not mcluded in the calibration 

We have discovered a new approach avoiding the above problems with 
20 correlations, and regression calculations, and being capable of being expanded 
automatically with use of a new product of ditferent quality 

The present invention provides a method of determining or predicting a value 
Px, of a property of a material X or a property of a product of a process from said 
material or yield of said process, which method comprises measuring the absorption 
25 D|X of said material at more than one wavelength in the region 600-2600nm, 

comparing signals indicative of said absorptions or mathematical fijnctions thereof with 
signals indicative of absorptions Djm at the same wavelengths or mathematical 
functions thereof for a number of standards S in a bank for which the said property or 
yield P is known, and choosing from the bank at least one and preferably at least 2 
30 standard with property Pm, said standard having the smallest average values of 
the absolute values of the difference at each wavelength i between the signal for the 
material and the signal for the standard to obtain value P^- w'^'"' averaging of 
said properties or yields Pm. when more than I standard is chosen 

The above method can be -performed without regression or correlation 
35 techniques, eg between the absorption at any wavelength of the material and the 
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property/yield 

Thus for the perf ormance of the method of the invention, a bank is prepared in 
which the NIR spectra are recorded at many wavelengths for a large number of 
standard materials, together with their propenies (or those of products obtained by 
5 processes therefrom) determined by alternative techniques e g. gas chromatography for 
chemical compositions and yields determined by known methods ©r viscosities by 
known mechanical methods The standards are chosen to cover the area in which the 
method is to be used, so for octane number determination, a range of gasolines is 
chosen of widely varying octane numbers, with different contents of lead, or other 

1 0 additives such as alkyi ethers and aromatics For the determination of properties of 

polyethylenes a range of polyethylenes is chosen of widely varying properties, e.g. with 
different contents of comonomer, or other properties such as molecular weight. For 
viscosity determinations for base oils, a range of base oils is chosen of widely varying 
viscosities. The number of wavelengths chosen may be 2-1000 e g 5-200 or 10-20 

1 5 such as 40-80 especially for oil refining/petrochemical operations as described below, 
or 5-100 or 10-80 such as 25-65 especially for use with processes which are 
polymerisation, oligomerisation or an organic reaction in which at least one of the 
reactant and a product is a functioaalised compound, or 10-80 such as 40-70 especially 
where material X is a composition comprising pan of a lubricating oil fraction from a 

20 distillation of oil. The number of standards can be at least 100 or 1000, or 100,000 up 
to 5 million depending on propeny(ies) cho.sen 

The wavelengths chosen may be at regular intervals such as each 1-50 or 10-50 
(especially for such processes like polymerisation oligomerisation and reaction as 
described above) or 1 5-35nm (or each l-5nm or each nanometre) or may be at irregular 

25 intervals e g with internals of l-200nin e g I- 100 or I -50 such as 2-50 or 4-50 or 10- 
60nm, which may be random or chosen because of a change in the shape of the spectral 
curve at that wavelength e.g. a peak, uough or shoulder or chosen by chemical or 
statistical criteria such as factor analysis. The wavelengths may be in the region 600- 
260Gnm. such as 800-2600nm, in particular 1500-2600 or 2000-2550nm, especially for 

30 oil refining/petrochemical operations as described below or 800-2600 eg 800-2000 
especially 1000-I800nm or 2000-2550nm for diene containing gasolines such as ones 
produced by cracking e g steam cracking The wavenumbers may be in the region 
16,600-3840cm-', e g 12.500 to j840cm-' in particular 6660-3840 or 5000-3900cm-l, 
or 12500 to 3840M 2500-5000 especially 10000-5500 or 5000-3900cm- 1 especially for 

35 oil refining/petrtfchemical operations as described below; corresponding frequencies in 
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Henz can be obtained by multiplying iliese wavenumbers by jxiO' OcnVsec. 
Wavelengths may also be in the region 600-2500nni, e g 900-2500nm such as 1000- 
2000nm, while the wavenumbers may be l6.600-4000cm-' such as I i 000-4000 or 
lOOOO-SOOOcm"'. both in particular tor the polymerisation, oligomerisation or organic 
5 reactions described above and beiou' Wavelengths may also be in the region 600- 
2600nm, e.g. 1000-2500nm but preferably 1500-2600 or 2000-2550nm, while the 
wavenumbers may be 16,600-3840cm-', e.g. 1 0000-4000cm- ^ e.g. 6660-3840cm-' or 
5000-3900cm"', especially for processes in which material X is a composition 
comprising part of a lubricating oil fraction from the distillation of oil, 
1 0 The signals eg absorptions (or derivatives) for the unknown sample are 

compared with the signals eg absorptions (or derivatives) at the same wavelength of 
the standards, and those standards chosen havmg the smallest differences The 
properties of those chosen standards are then averaged to determine the propeny of 
the unknown sample The absorptions at more than one wavelength may be chosen. 
1 5 e.g. 2-1000 such as 5-100 or 10-20 Other methods of signal processing apart from 
derivatives such as Fourier transformation may be used in a similar way 

In the method of the invention the standards chosen are those with the smallest 
average values of the absolute ditTei ence at each wavelength i between the signal 
exemplified by absorption/optical density (or a derivative thereof) Dj^ for the unknown 
20 material and the corresponding signal eg absorption/optical density (or derivative 

thereof) Djm for the standard The averages may be in respect of the mean value of 
Djx-Djm (whatever its sign i.e absolute difference), or (Dix'^jm)^ and may be the 
simple mean value or the differences may be weighted to take accoimt of the different 
sensitivity of the absorption to the propeny at that wavelength or the different 
25 sensitivity of the spectrometer at tliat wavelength For each standard in the bank of 
standards for the type of material in question, the average difference is found as 
described and the standard or standards with the smallest average differences chosen, 
e.g. at least I but preferably at least 2 such as upto 1000 smallest such as I (or 2)- 1 00 
or 1 (or 2)-20 but is particular I (or 2)- 10 and especially 2-6 smallest Advantageously 
30 the average differences chosen and hence the standard (or standards) chosen for 
the property or yield wanted are such that in relation to the unknown material X and 
each chosen standard the following functions is met 
'xm experimental error 



35 
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wherein "S the proximity index and is defined by |2(xin) = I{D,x-Dim)2 and the 
experimental error is in determining said property or yield in the standard The value 
Px of the propeny or yield is the same as propeny or yield or the average Pm if 
more than one standard Sm is chosen 
5 In order to aid the choice of the appropriate standards, especially in relation to 

a large number of wavelengths for a complex unknown mixture, it js preferred to limit 
the choice to those defined by means of a minimal index. For the chosen standard the 
minimal index is at least the same as the differences between the absorptions of the 
unknown and the standards. Mathematically, this may be expressed as i^abSi^M 

1 0 where iM is the minimal index for the propeny, and iab is a measure of the deviation 
(called the proximity index) at all the chosen wavelengths between absorption of the 
unknown and a chosen standard b That measure is defined by 

i(ab)2 = li (Di,-Dib)2 (I) 
where Dja is the optical density (or absoi bence) of unknown a at wavelength i (or a 

1 5 derivative thereof e.g. a first, second or third derivative of that density), and Djb is the 
optical density (or absorbence) of standard b at that wavelength i (or a derivative 
thereof e.g. a first, second or third derivative of that density). The value of D) is the 
optical density or the optical density dilference with respect to the baseline of the 
spectrum at that wavelength, or the baseline interpolated between 2 wavelengths on 

20 either side thereof If desired signals corresponding to other mathematical functions of 
the absorption eg after Fourier transformation or spectral subtraction or division may 
be used to provide corresponding proximity and Minimal Indices. 

If desired instead of the optical density Dj a normalized density Wj may be used 
where Wj =Di/lDi. This normalization avoids errors due to small electronic 

25 fluctuations in the apparatus and compensates for small differences in the optical path 
between the optical cells.. In this case the proximity index is defined by 
I(ab)2 = Ij (Wia-Wib)2 (2) 
The indices can be weighted as desired for increasing resolution One approach 
is to define the indices as follows. 

30 I(ab)m = XAbs value (Xij,-Xjb)'^/ajn (3) 

where Xj is Dj or Wj or a mathematical combination thereof, o\ is the standard 
deviation of X for the set of samples considered (at that wavelength) and each of m 
and n which are the same or diflerent is weighting factor which is positive but can be a 
<whole number ora fraction. Other variants can be used with other weighting factors 

35 such as those involving the spectral experimental error e;, where ej is the 
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reproducibility of the spectral measurement at wavelength i The choice between the 
different options for the weighted indices may be dictated by numerical efficiency. 

The reproducibility of the experimental measurements in the standards may be 
at least 90% or 94% or 95%. The minimal index may be obtained from a reference 
5 standard samples set according to the following procedure, hereafter called the 

Minimal Index Procedure The NIR spectra for 2 standard samples A and B and their 
property P e.g. Octane Number, or viscosity e g for a polymer, or density e g for a 
lubricating oil fraction are determined By means of equation ( I ), (2) or (3), the value 
of the proximity index i^b is determined via the absorptions at a series of wavelengths; 

1 0 this index is applicable to the difference in properties Pa-Pb called EPj^b 

This process is repeated with other pairs of standards c and d, e and f etc to 
obtain a series of Proximity Indices ic^^ etc with corresponding propeny differences 
EP(.(j etc For different values of a parameter L which is greater than the indices i^b 
etc. the corresponding values of EPj^b etc are averaged to give an average EPy for that 

1 5 value of L, the different values of EPjj+ta/VK are then plotted on a graph against L o 
is the accuracy of the propeny determination and K is the number of pairs of samples 
for which i^b is inferior to a given L t is the Student factor at a given level of 
confidence. The intercept is then measured between the curve obtained and a line 
usually horizontal which is the reproducibility of the propeny level at an appropriate 

20 confidence interval e g 90% or more usually 95%; the abcissa ponion of the intercept 
gives the minimal index 1,^;^, which is the minimum value of i^b for which Pa=Pb 
within the frame of experimental error 

From this minimal index by Procedure I. the standards can be chosen which 
have values of i-j^b ^ i-min where m this ca.se a is the unknown and b is a standard, as 

25 in this case the difference between Propeny a and Propeny b is less than or equal to oV 
2. where a is the experimental error in measuring the propeny Then from the 
propeny P value or values of the chosen standard, the propeny of the unknown is 
obtained directly or by averaging those values, usually the arithmetic mean, but 
optionally with weighting. 

30 The method of the invention may be used to determine more than one Property 

P at once, e.g. at least 2, such as 1-30 e g 2-10 properties at once. Each propeny of 
the standards has a particular unweighted, minimal index, which may lie in the region 
0-10-10 e.g. 10-1 to 10-'^ or 10-2 to 10-^, in particular lO''^ (or 5x10-7) to 5x10-^ for 
Minimal Indices derived from absorbancies; corresponding Minimal Indices may be 

35 obtained for othty signals/functions If the Minimal Index chosen is the smallest for all 
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the properties desired, then ihe same one may be used for all the properties and the 
standards chosen will be suitable for all the properties The Minimal Index for each 
property may be used separately, with difterent numbers of standards chosen for each 
property (assuming different Minimal Indices) If desired the same Minimal Index may 
5 be used, which is not the smallest, resulting in some of the chosen standards (with a 
higher Minimal Index) giving some properties of high accuracy and some (with a lower 
Minimal Index) giving some properties of less high accuracy. 

The property to be determined may be of the sample being analyzed or a 
product obtained from that sample e g. a product of blending, cracking, separating or 
1 0 polymerising the sample, as the propeny value obtained is derived from the standards, 
and they will have been determined as needed for the eventual use Our EP304232 and 
305090 referred to above describes such techniques when applied to use of NIR with 
correlation to blending, separatini! or cracking operation; the same principles apply in 
the present method. 

1 5 If the density of the standards in the data bank is sufficient to have \^a.h<i^ min 

as is usually the case, the above procedure is vei^ satisfactory. But there are occasions 
when the bank is incomplete, because of shortage of data of propenies in a particular 
area i.e. a low density of standards or the sensitivity of the propeny to changes in 
absorption is so small, that a very small Minimal Index is required and there may be 

20 few standards with proximity indices meeting it It is possible simply choose a larger 
Minimal Index with e.g. \-5 times such as 1.5-2 times the Minimal Index; the results 
may be less accurate than those from a smaller minimal index. 

However, a more accurate approach with a low density of standards involves a 
special densification process of Procedure 2, in which random or semi random 

25 densification of the neighbourhood of the unknown is achieved by generation of 
synthetic standards, based on standards already in the bank Each new synthetic 
standard may be obtained from combinations of standards taken at random from the 
bank but preferably it is obtained from the other standards by the constraint of 
choosing only a mixture of N standards for which 

30 (Min)Cj - uj <C\j <(Max)Cj+uj (4) 

and ICij = 1 (5) 
where Cjj is the fraction of component j in the sample \ 

Min Cj is the minimum amount of j in the initial calibration mixture i e. 
'Standards in the bank or in the samples for which the method is to be used, and 

35 Max C; \t the maximum amount of; in the initial calibration mixture i.e 
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Standards in the bank or in the samples for which the method is to be used, and 

uj is usually between 1 and 0.01 or I and 0 05 preferably between 0.5 and 0. 1 
and can be fixed for each property. 

The constraints over the choice of such mixtures of N standards can also be 
5 equally fixed in the spectral area from which the samples will be drawn in order to 
remain in the areas of similar chemical nature 

The number of samples effectively drawn into the bank in this densification can 
be of several thousand generally 1000-2000. The calculation time is extended without 
significant deterioration in the results, If no further neighbours are found, the trawl of 
1 0 new samples drawn in is enlarged 

The spectrum of each mixture Is calculated by the combination of the spectra of 
the standards used according to the formula 



where Pj is the propeny of component j 

In the case of non linear additive propemes. appropriate mixing factors can be 
20 applied e g by blending factors or similar for density and viscosity 

Having obtained the spectrum and the properties of the synthetic mixtures, 
these can be used as "standards" to help determine the properties of an unknown 
sample in the same way as a conventional standard 



25 Procedure 3 may be used as follows The Q nearest samples to unknown X can be 
found from a selection from the bank samples for which the proximity index to the 
unknown sample is (V) X i^^in) where v is 0 I < v < 10, (8) preferably 0 5<v<2 or 1< 
V <5. Then by the method of least squares is found a generally linear combination of 
the standard products, which are the Q nearest samples, to reproduce the spectrum of 

30 X according to the equation 



where Cr is the coetficient for sample R in the total Q and Sr is the spectrum of 
sample R. The coefficient Cr which can be normalized to Cr^ = I or not and/or 
optimized by the least squares route, allows an estimation of the property Px according 



15 




Instead of using either of the two above approaches, 1-7, a third type 



Sx= ICrXS, 



(9) 



35 



to the equation 



8 
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Px= iCRXPf^ (10) 
where Pr is the propeny of sample R 

The eventual size of the estimation error can be derived by application of 
Gaussian theory, also called the propagation error (see Eq 10) 
5 The above third approach can only be applied if the product X is situated inside 

the maximum extension of the standard products defmed by equatjon (8) i.e v^ithin the 
range of bank samples defined in equation (8). If this is not the case, X is outside the 
field of the actual bank of products and escapes from the area of knowledge of the 
method into the area of learning 

1 0 The densitkation process described in relation to equations 4-7, or 9 or 10 is 

usually applied to the method of the invention involving no correlation or regression 
techniques. However, if desired the densification process may be applied to increase 
the number of "standards" for consideration in an NIR analytical technique involving 
the correlation on regression techniques as described above e g. MLR The present 

1 5 invention also provides a method for adding an extra synthetic standard to a bank of 
known standards, each of which relates at least one absorption in the 600-2600nm 
region (or a signal indicative thereof or of a mathematical flinction of said absorption 
eg a derivative thereoQ of a known material to a known properly related to that 
material, which method comprises choosing from the bank at least 2 standards for 

20 which equations 4 and .S above are met. considering mixing the chosen standards in at 
least one proponion to produce at least one mixture for use as a synthetic standard, 
and estimating the spectrum and property/yield of said mixture according to equation 6 
and 7 respectively 

The spectrum and propeny/yiekl of each "mixture" can then be added to the 
25 bank and used to develop models through the known correlation/regression approach, 
e.g. as described in the above mentioned patents. 

The method of the invention may be applied from the spectrum of a material to 
determine at least one physical, chemical, physicochemical and/or rheological property 
of that material, which may be a product of a chemical or physical or separation 
30 process, or which may be a feed to such a process, or the method can be used to 
determine at least one of said properties of a product of that process from the 
spectrum of at least one feed to that process, or to determine the yield of at least one 
product of that process. Each of the feed (or feeds) or products to the process may be 
t a solid liquid or gas preferably at least one feed or product is a liquid. 
35 Thus the method may be used tor the physicochemical determination or 
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prediction in relation to at least one teed or product used in or obtained by an 
industrial process of the refining of oil and/or in petrochemical operations The 
process may be a hydrocarbon conversion or separation process, preferably a 
reforming or catalytic cracking or hydrotreatment process or distillation or blending. 
5 In particular it may be used for determination of at least one property of a feed and/or 
the prediction and/or determination of at least one propeny and/or,yieId of product 
from a number of different processes such as processes for separating petroleum 
products such as atmospheric distillation vacuum distillation or separation by 
distillation, under pressure greater than atmospheric, as well as thermal or catalytic 

10 conversion, with or without panial or total iivdrogenation, of a petroleum product, 
such as catalytic cracking e.g tUiid caialyiic cracking (FCC), hydrocracking, 
reforming, isomerization. selective liydrogenation. viscoreduction or alkylation 

Of panicuiar value is rhe use ot the method in blending operations involving the 
prediction and/or determination ot'ai least one propeny of a blend of liquid 

1 5 hydrocarbons (optionally with other atlditives such as alkyi ethers), this method 

including or not the determination for each constituent of the blend of a blend index for 
the property considered. In this method as applied to blending, the blend indices can 
be obtained simply by calculation and without the need for preparation of standard 
physical mixtures other than those contained in the databank. The blend indices can be 

20 combined linearly or non linearly within the fields of stability to determine from the 

value of this combination a value for at least one property of the blend obtained The 
blend may be made by mixing at least 2 of butane, hydrogenated steamcracked 
gasoline, isomerate, reformate. MTBE di TAME. FCC derived gasoline This process 
may be repeated with numerical addition of other constituents separately to the liquid 

25 hydrocarbon base to determine a series of blending indices and then determination 
from these indices of the properties uf the multi constituent blend (see eg. Ex. 2 
hereafter). 

Examples of properties that can be determined and/or predicted include the 
following: for automobile fuels/gasolines, at least one of the Research Octane Number 
30 (RON). Motor Octane Number (MON) and/or their arithmetic mean, with or without 
lead additive and/or the methyl ten, butyl ether or methyl isoamyl ether and/or benzene 
content: 

For automobile fuels/gasolines, at least one of the vapour pressure, density, 
volatility, distillation cur\'e, e g percentage distilled at TO^C and/or IOO°C, oxygen 
35 content or benzene or sulphur content, chemical composition and/or gum content e g 
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expressed in mu/lOOml. and/or susceptibility to lead (these propertres are particularly 
determined for use in blending operations) 

For diesel fuels or gas oils, at least one of the cetane number (e g motor 
measured), cetane index, cloud point, "discharge point", filterability, distillation curve, 
5 density e.g. at 1 5°C, flash point, viscosity e.g. at 40°C. chemical composition, 
sensitivity to additives and percentage of sulphur; 

For distillation products from cinde oil e.g. under atmospheric pressure at least 
one of the density, percentage of sulphur, viscosity at 100°C, distillation curve, paraffm 
content, residual carbon content or Conradson carbon content, naphtha content, flash 

10 point for petrol, cloud point for gas oil e.g. light gas oil and/or viscosity at 100°C 

and/or sulphur content for atmospheric residues, and yield for at least one of the cuts, 
gasoline (bp 38-9.S°C), benzine (bp 9.S-I49°C) naphtha bp 149-1 75°C. jet fuel bp 175- 
232°C, light gas oil bp 232-j42°C, heavy gas oil bp ?42-?69°C. and atmospheric 
residue greater than 369°C 

1 5 For at least one of a feed or a product of a process of a catalytic cracking e.g. 

FCC process, at least one of the density, percentage of sulphur, aniline point, gas oil 
index, gasoline index, viscosity at IOO"C. refractive index at 20°C and/or GCC. 
molecular weight, distillation temperature e g 5>0% distillation temperature, percentage 
of aromatic carbon, content of total nitrogen and factors characterizing the suitability 

20 of the feed for the cracking e g KUOP. crackability factor, cokability factor, and yield 
e.g. of gas. gasoline, gas oil or residue Thus there may be determined the yields 
and/or properties of the difterent products obtained by distillation of the cracked 
products, such as RON and/or MON, clear or leaded for the gasoline cut and the 
viscosity at 100°C for the distillation residue. 

25 For at least one of a product or a feed of a catalytic reforming process, at least 

one of the density, distillation temperature and/or chemical composition (expressed as 
a percentage) of saturated linear hydrocarbon. isoparatVins. naphthenes. aromatics and 
olefins. 

For at least one of a product or a feed of a process of hydrogenating gasoline at 
30 least one of the density, distillation temperature, RON and/or MON. clear or leaded 
vapour pressure, volatility, chemical composition (expressed as a percentage) of 
saturated linear hydrocarbons, isoparatlins, naphthenes, aromatics e g benzene, and 
mono/di substituted benzenes, olefins e g. cyclic and non cyclic olefins, diolefins. the 
nialeic anhydride index, and yield eg. of at least one of the products obtained. 
35 The method of the invention may also be used with chemical reactions in which 
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ai least one product is a hydrocarbon, and none of the feeds or products contains an 
element other than carbon or hydrogen The hydrocarbon which may be gaseous or 
hquid at 25°C. Such reactions may involve as feed or product at least one olefin or 
acetylene e.g. linear or branched, aliphatic or cycloaliphatic olefm with an internal or 
5 external ethylenic unsaturation. preferably of 2-20 carbons especially 2-8 carbons for 
alkenes or alkynes (such as ethylene, propylene, butene I or 2, isobiitene, isopentene) 
or acetylene, and 5-8 carbons for cycloalkenes e g cyclohexene. The feed or product 
may also be an aromatic hydrocarbon e g benzene or naphthalene, optionally 
substituted by at least one (e g. 1-3) alkyi or alkenyl group e.g. of 1-20 carbons, such 

10 as 1-6 carbons, especially methyl, ethyl or isopropyh examples are benzene, toluene 
xylene, cumene and styrene The feed or product may also be a non aromatic 
hydrocarbon, e g linear or branched aliphatic or cycloaliphatic with e g 1-20 or 5-8 
carbons respectively, preferably I -6 carbons and 6 or 7 carbons respectively, examples 
are methane, ethane, propane, n-butane. isobutane, and cyclohexane. The feed or 

1 5 product may also be a diene. conjugated or unconjugated, aliphatic or cycloaliphatic 
with e.g. 4-20 carbons or 6-20 carbons respectively, examples are butadiene and 
isoprene and cyclohexadiene. Examples of the reactions are hydrogenation {e.g 
butadiene to butene- 1 or 2 or cyclohexene to cyclohexane) dehydrogenation (e.g. 
ethane to ethylene or ethyl benzene to styrene). isomerisation (e g butene- 1 or -2 to 

20 isobutene, or pentene-1 to isopentene) aikylation (e g benzene with ethylene to form 
ethylbenzene and/or styrene. or isobutene with butane to form iso octane), and 
cracking 

In addition to the use in petrochemical operations, the method is of wider 
application and may be applied in the pharmaceutical industry such as the production 

25 of pharmaceutically active compounds for use as medicines e g by fermentation, and in 
the perfumery industry for making perfumes and fragances. especially in their blending 
and control thereof The method may also be used in the food industry e.g. in brewing 
to control fermentaion processes, in fermentation to make wine and quality control 
thereof, and control of food production e.g. sugar and water content in fruit juice and 

30 in control of maturing processes for fruits and vegetables. In each case the method 
may be applied to determine a propeny of the sample tested or product from that 
sample e.g. a fermentation or blended product preferably on line and especially with 
continuous feed back from the results to control the production process 

The knowji correlative techniques for modelling physicochemical properties 

35 based on NIR sjiectra have di.sadvantages, the chief of which is the need for 
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establishing a strong correlation between the spectmni and the propeny, and their 
difficulty in dealing with positive or negative synergy between components 
contributing to that property. In the case of high density polyethylene one multi linear 
regression model in respect of density gives a coefficient of correlation that can on 
5 occasion be insufficiently high, so as to give problems in a polymerization process 
based on it. 

The method of the invention can avoid the above problems with correlations, 
and regression calculations, and can be capable of being expanded automatically with 
use of a new product of different quality 

"•0 The present invention also provides a method of determining or predicting a 

value Px. of a property of a material X which is a feed to a process or a property of a 
product of a process from said material or yield of said process which method 
comprises measuring the absorption (D,m) of said material at more than one 
wavelength in the region 600-2600nm. comparing signals indicative of said absorptions 

15 or mathematical functions thereof with signals indicative of absorptions Djm at the 
same wavelengths or mathematical functions thereof for a number of standards S in a 
bank for which the said property or yield P is known, and choosing from the bank at 
least one, and preferably at least 2 standard Sm with properties or yield P^ said 
standard Sm having the smallest average value of the absolute difference at each 

20 wavelength i between the signal for the material and the signal for the standard S^, to 
obtain value P.^^ with averaging of said properties or yield Pm. when more than 1 
standard 5,^, is chosen, and wherein said process is at least one of a polymerization, an 
oligomerization or an organic reaction in which at least one of the reactant and a 
product is a funciionalized compound The above method can be performed without 

25 regression or correlation techniques, e g between the absorption at any wavelength of 
the material and the property/yield 

The present invention also provides a method for adding an extra synthetic 
standard to a bank of known standards, each of which relates at least one absorption in 
the 600-2600nm region (or a signal indicative thereof or of a mathematical function of 

30 said absorption eg a derivative thereof) of a known material to a known property 

related to that material, wherein .said property is of said material, which is a feed to a 
process, or product of said process, or yield of said process, which method comprises 
choosing from the bank at least 2 standards for which equations 4 and 5 above are met, 
considering mixing the chosen standards in at least one proportion to produce at least 

35 one mixture for ase as a synthetic standard, and estimating the spectrum and propeny 

13 
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of said mixiui e according lo eiiuation o and 7 i cspeciively. said process being at least 
one of a polymerization an oliuomerizaiion or an organic reaction in which at least one 
of the reactant and a product is a funciionalized compound. The spectrum and 
property of each "mixture" can then be added to the bank and may be used to develop 
5 models through the known correlation/regression approach, eg based on that in the 
patents described above. 

The method of the present invention is applicable to chemical reactions, which 
may be polymerisations or oligomerisations. or alternatively reactions in which at least 
one of a reactant and a product is a functionalised compound. In the chemical 

1 0 reactions each of the feeds and the products may be a solid, liquid or gas, preferably all 
the feeds are liquids and/or gases, and preferably all the products are liquids and/or 
solids, especially liquids 

Examples of polymerisatioii.s are condensation and addition polymerisation 
Condensation polymerisations may pioduce thermosei polymers, such as phenolic 

1 5 novolac or resole resins curing with oi without curing agents like hexamine. or 

polyurethanes. or thermoplastic polyiueis such as polyamides. e.g. polylactams such as 
Nylon-O and polymers from polyaniiiit.s and polycarboxylic acids e g poly 
hexamethylene adipate, and polyesitrs. such as those from diols e.g. aliphatic diols and 
organo di carboxylic acids e.g. aromatic or aiyl bis (alkylene) dicarboxylic acids, such 

20 as poly ethylene terephthalaie Addition polymerisations tend to produce 

thermoplastic polymers, and may be thermal or free radical or catalysed reactions e g 
with Bronsted or proton acids or metals, especially transition metals. Examples of 
such polymerisations are those involving polymerisation at an olefmic double bond or 
ring opening of an epoxide or episulphide Tlie olefmic double bond is preferably a 

25 vinyl group CH2=C- and may be in a iivdrocarbon e g an olefin especially an alkene 
such as one of 2-12 carbons especially ethylene alone or mixed with at least one alpha 
olefin of 3-12 carbons (especially in amount ot'O 5-30% by weight based on total 
olefins) such as propylene, butene- 1 , 4-inethYl-pentene- 1 , hexene- i , octene- 1 or 
styrene; copolymers of such olefinic hydrocarbons, especially ethylene, with non 

30 hydrocarbon comonomers e g. esters with olefinic groups such as vinyl esters e.g. vinyl 
acetate or alkyl(meth)acrylate or vinyl chloride may also be made Addition 
polymerisation of iso olefins e g. of 4-b carbons such as isobutene alone or with other 
comonomers such as butadiene is included, as in addition polymerisation of olefinic 
non hydrocarbon monomers such as vinyl esters e g of 3-20 carbons especially 4-10 

35 tarbons such as vinyl acetate and propionate, and alkyl(meth)acrylates wherein the 
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alkyi group has 1-20 carbons, especi<illy 1-4 carbons for solid polymers e g polymethyl 
methacrylate, and 4-20 carbons for polymers Ibi use as pour point depressants and VI 
improvers e.g. poiydodecyl acrylaie and nietliaci-ylate and copolymers with 2-10 
monomers of different alkyl chain lengths Vinyl chloride homopolymers and 
5 copolymers e g with vinylidene chloride may also be made 

The method may also be used for ring opening reactions such as reactions of 
epoxides, episulphides or cyclic immes with organic compounds containing at least one 
active hydrogen such as compounds with at least one OH, NH or SH group, such as 
alcohols, phenols, primary or secondai y amines or thiols Alcohols e g. of 1-30 

1 0 carbons such as 2-6 carbons (e.g. butanul) especially alkanols or cycioalkanols are 

preferred. The epoxide is usually carbons e.g. ethylene oxide, propylene oxide, 

butylene oxide or cyclohexane oxide, while the episulphide and cyclic imines are 
preferably the corresponding analogues e g. ethylene imine and ethylene sulphide. 
In the case of polymerisation the method may be used to estimate the 

1 5 properties of the polymer made from the NIR spectrum of the feedstock (under 

constant conditions) oi from the NIK spcctnim of the product Examples of properties 
are number and weight average molecular weights, and the molecular weight 
distribution, viscosity e.g. at lOU^C", tluidiiy aidex, density, and chemical composition 
e.g. percentage of at least one monomer or comonomer in the polymer percentage of 

20 unsaturation e g ethylenic type, or side chain grouping, e g methyl, crystallinity, 

rigidity, tlow parameters, draw strength at the tlow threshold, free cracking resistance 
and shock resistance. In addition for polyisobutenes. the property may also be content 
of butene-l, and light and heavy polyisobutenes and unsaturation expressed in groups 
per litre and maleinisation index (or succinylation ratio) (sensitivity to Diels AJder 

25 reactions) as well as particular types of unsaturation e g vinylidene CH2=C- VfN, tri 
(CH3-C(CH3)=CH<) TRIl, tii 2 (TRI2) (CH3-CH=C-) TRITOT (R-CH=C<), 
TETRA (>C=C<). For polyolefms e g polyethylene, other properties include 
percentage of comonomer, volatile compounds and degree of conversion. For 
polyalkylenoxylated compounds e.g. ethylene oxide condensates e g with alcohols, the 

30 method may be used to monitor the degree of conversion or the amount of alkylene 
oxide consumed, as well as the quality of the product e g content of groups derived 
from at least one epoxide or the di.stribuiion of those groups in the polymer chain, the 
product weight and number average molecular weight and its distribution, proportions 
'of low and high molecular weight products (e g 1 50-600 or 600- 1 .SOOO such as 5000- 

35 12000 respectively) Hydroxyl index (or mean number of hydioxyl groups per 
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molecule), percentage ot" primary secondary' and tertiary iiydroxyl groups, allylic or 
propyienic type unsaturation. or impurity content 

The method is of especially value in the polymerisation of ethylene alone or 
with at least one alpha olefin as described above The process is usually catalysed by 
5 at least one transition metal catalyst especially of Group IVA, VA or VIA, of the 

Periodic Table, such as litaniuni, zirconium, vanadium and/or chromium The catalysts 
may be organometallic (including D complexes), especially with the above transition 
metals, and may be in the presence of at least one organo aluminium cocatalyst as in 
Zeigler Natta catalysts Non organometallic catalysts such as chromium oxide may be 
1 0 used. The catalyst may be unsupported or supported e.g. on silica and/or alumina. 

The method may also be applied to organic chemical processes, which are not 
polymerisations (including oligomerisation); thus processes involving only monomeric 
starting materials and products are suitable In particular these include processes in 
which at least one of a reactant and a product is a functionalised compound i.e. is not a 
1 5 hydrocarbon but contains at least one functional group, e g with at least one atom 
other than carbon and hydrogen, in panicular at least one oxygen, nitrogen, sulphur, 
phosphoais. or halogen e g chlorine, bromine, iodine or tluorine atom, especially 1-3 
such atoms in the compound. The functional group may be an alcohol, phenol, thiol, 
primary secondary or tertiary amine, aldehyde, ketone, ester, acid, amide, nitrile or 
20 ether or sulphide, or aromatic or aliphatic halide 

In particular the process may be a hydration such as an olefm to an alcohol 
(e g. ethylene or propylene to ethanol or isopropanol respectively) dehydration such as 
an alcohol to an olefin (e g tert butanol to isobutene) etheritlcation such as reaction of 
an alcohol or phenol with an olefm (e g leit butanol wiih isobutene to form Methyl 
25 tea butyl ether) or reaction ot an oleliii with water (e g ethylene to diethyl ether), 
esterification such as reaction of a cnrboxylic acid (or derivative thereof e g. acid 
chloride) with an alcohol e g alkanol of 1-20 carbons) or with an olefin (e g. ethylene, 
propylene or n or isobutene), such as reaction of acetic acid with ethylene to form ethyl 
acetate or with dehydrogenation) vinyl acetate The process may also be an oxidation 
30 e g. an alcohol or aldehyde to an acid such as methanol to formic acid, or a 

hydrocarbon to an alcohol or ketone or an acid e g naphtha to acetic acid or methane 
to formic acid or cumene to acetone and phenol, an ammoxidation e.g. an aliphatic 
substituted olefm (with optionally ?>-(> carbons) to a nitrile such as propylene to 
acrylonitrile, or a cai bonyiation of an oletin or an alcohol to form a carboxylic acid 
35 and/or anhydride^ such as the reaction of methanol with carbon monoxide to form 
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acetic acid and/or anhydride. 

In a further embodiment the pi esent invention provides a method of 
determining or predicting a value Px. ot a property of a material or a property of a 
product of a process from said material or yield of said process which method 
5 comprises measuring the absorption (Djx)of said material at more than one wavelength 
in the region 600-2600nm. comparing signals indicative of said absorptions or 
mathematical functions thereof with signals indicative of absorptions Djm at the same 
wavelengths or mathematical functions thereof for a number of standards S in a bank 
for which the said propeny or yield P is known, and choosing from the bank at least 

10 one. preferably at least 2 standard Sm with property Pm having the smallest average 
value of the absolute difference at each wavelength i between the signal for the 
material and the signal for the standard S,,, to obtain P,-,^ with averaging of said 
properties or yields Pm when more than I stnndaixi Spi is chosen and wherein said 
material X is a composition comprising part of a Uibricaiing oil fraction obtainable 

1 5 from a distillation of oil e.g. a vacuum distillation of oil 

The above method can be performed without regression or correlation 
techniques, e.g. between the absorption at anv wavelength of the material and the 
property/yield 

The densification process described in relation to equations 4-7, or 9 or 10 is 

20 usually applied to the method of the invention involving no correlation or regression 
techniques However, if desired the densification process may be applied to increase 
the number of "standards" for consideration in an NIR analytical technique involving 
the correlation on regression techniques as described above e.g. MLR. The present 
invention also provides a method for adding an extra synthetic standard to a bank of 

25 known standards, each of which relates at least one absorption in the 600-2600nm 

region (or a signal indicative thereof or ota mathematical function of said absorption 
eg a derivative thereoO of a known material to a known property related to that 
material, wherein said propeny is of said material or is of a product or yield of a 
process from said material, which method comprises choosing from the bank at least 2 

30 standards for which equations 4 and f> above are met, considering mixing the chosen 
standards in at least one proportion to produce at least one mixture for use as a 
synthetic standard, and estimating the spectrum and property of said mixture according 
to equation 6 and 7 respectively and wherein said material X is a composition 
comprising part of a Jubricating oil fraction obtainable from a vacuum distillation of oil 

35 e.g.a vacuum distillation of oil 
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The spectaim and property ol each "mixture" can then be added to the bank 
and may be used to develop the models tlirouyh known correlation/reuression 
approach, e.g. as described in the above mentioned patents. 

This embodiment of the method of the present invention is applicable to various 
5 petroleum hydrocarbon tractions, which comprise pan (and only pan) of a lubricating 
oil fraction e.g. from a vacuum distillation of oil after removal of materials boiling 
above 310°C (under atmospheric pressure). Such fractions include the partly purified 
lube cut from the distillation, e g. after at least one of the steps of dewaxing and 
dearomatizing and preferably both, (as in lube base oil) and the partly purified vacuum 
10 distillation residue e g. after at least one of the steps of deasphalting, dewaxing and 

dearomatizing, and preferably all 3 (as in bright stock). Such fractions also include the 
aromatic extract of the lube oil cut or distillation residue, or a wax separated 
therefrom. 

The method is preferably applied to lube base oils or bright stock The base oil 

1 5 may be a 100-600 neutral or solvent or BS oil e g. 100. I 50. 200. 300. 400 or 500 
neutral oil or BS solvent It may have at least one of and preferably all of the 
following properties a density at 1 5"C ot O SO-U 95kg/l e.g. 0 85-0/92kg/l. a kinematic 
viscosity at 40°C of 10-1 OOOcSi e u 1 .^-70UcSi. and at 1 00°C of 0 5-50cSt e.g. 1- 
40cSt. a Flash Point of 180°C mm e u l'^)0"C min, a pour point of 0°C maximum e.g. - 

20 5"'C or -7°C maximum and a Viscosity Index of SOmin e g 90min The base oil may 
be present alone, or may be mixed with the aromatic extract as in process oils, which 
may have at least one of, and preferably all of the following properties, a density at 
I5°C of 0.95-l.lOkg/l. eg 0 97-1 OGkg/1. a Kinematic Viscosity at 40°C of at least 
30cSt e.g. at least 37cSt. and at 100°C of at most 50cSt e g. at most 45cSt and a Flash 

25 Point of at least 185°C e g I90°C mm The base oil may also be present mixed with at 
least one wax e g. in amount of 0-50% such as i-40% or 1 5-35% by weight as in 
"slack wax", the mixture of oil and solid wax separated in the dewaxing step, or waxes 
as in the residue from the dearomatization step 

The base oil may also be mixed with at least one non hydrocarbon additive to 

30 boost its effectiveness for lubricant use Types of additives which may each be present 
in amounts of 0 01-10% by weight (based on the weight of base oil) e g 0 1-1% are (i) 
detergents/dispersanis such a.'; alkyi phenates and/or alkyl aryl sulphonates (ii) 
antioxidants such as phenol derivatives, (iii) viscosity index improvers and pour point 
^ depressants, such as alkyl poly(meih)<ici ylate homo and especially copolymers, styrene 

35 butadiene polymers and polyisobutvlene (iv) ami corrosives, such as sulphur 



wo 96/11399 



PCT/GB95/02353 



compounds, zinc sulphophosphates and dithiophosphates. and (v) solid or liquid 
lubricity additives, such as graphite, molybdenum disuiphide and silicones. 

The method may also be applied to the aromatic extract resulting from the 
extraction of aromatics (e.g. with furfural) from the lube cut of the vacuum distillate or 
5 the deasphalted vacuum residue. This aromatic extract is different from the base oil as 
it contains a much higher amount of aromatics. such as benzene, toluene and xylenes, 
and higher molecular weight aromatics e.g. of at least 30 carbons than the base oil. 
The aromatic extract may be used alone or mixed with an amount of base oil to form 
process oil. 

0 The method may also be applied to solid or liquid paraffins or waxes e.g. as 

separated in a dewaxing step from the lube cut or the deasphalted residue. The wax 
may be mixed with base oil as in slack wax. or substantially free of base oil and may 
then if desired be further purified to produce a paraffin. Waxes may be used 
industrially while paraffins may be used for food and cosmetic uses 

1 5 The method is preferably applied for process control in a part of a refinery 

producing lubricants and by products therefrom, but may also be used for identification 
of unknowns e.g. for "finger printing" oils such as formulated oils. 

Examples of properties that can be determined/estimated for the various 
materials are as follows Where the material is a base oil (or formulated oil) the 

20 property may be at least one of the density, sulphur content. Flush Point. Flow Point, 
kinematic viscosity at 40°C and at IOO"C, Viscosity Index, aromatic carbon content, 
Polycyclic Aromatic hydrocarbon content, nitrogen base content, and inflammability 
according to Pensky Maiiens °C When the material is a crude paraffin or slack wax, 
the property may be at least one of the density, viscosity e g at 40°C or 100°C and oil 

25 content. When the material is a process oil, the property may be at least one of the 
density, sulphur content. Polycyclic Aromatic hydrocarbon content, viscosity e.g. at 
40°C or IOO°C and the Flash Point e.g. Cleveland Flash Point 

In each of the above processes the property or yield of a product determined or 
predicted by the method of the invention can be compared to the desired figure and 

30 notice taken of any deviations by adjusting the parameters of the process e.g. flow 
rates proportion or nature of feed(s) (e.g. via operation of control valves) and/or 
temperature/pressure etc to bring the property or yield back to the desired figure This 
control of the process, which may be a blending, separation or chemical e g 
polymerisation process, is usually performed with a micro computer which is linked to 

35 the spectrometer and also performs the search for the standards Sm The inline control 
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of the process is very efficient and very fast 

The present invention also provides an apparatus suitable for carrying out the 
method of the invention comprising an infra red spectrometer and a computer wherein 
the infra red spectrometer is linked to the computer programmed in such manner that 
5 the propeny or yield may be determined continuously and in real time. The 

spectrometer is suitable for measuring spectra in at least partly in the 600-2600nm 
wavelength range and can be linked to a signal processing device to allow numerical 
treatment of the spectrum, preferably by Fourier Transformation. The spectrometer 
receives at least one signal from a vessel containing product or from a feed or product 

10 line. The information obtained can be used as an information vector for the computer 
which is programmed to determine the property or yield e g via calculations on the 
proximity indices in relation to siandard.s Convenienily in relation to a process, the 
computer may be used in a closed loop feed back or feed forward control system for 
controlling processing equipment e u ciiangmg the process parameters in response to 

1 5 variations in the propeny and/or yield of product from the desired value, from 

measurement of more than one absorptions in the NIR spectrum of the product and/or 
feed 

The present invention also provides a computer programmed to perform the 
method of the invention of determining or predicting the value Px of the propeny or 

20 yield, or programmed to perform the method of the invention of adding an extra 

synthetic standard to the bank of known standards. The apparatus for use with the 
former method of the invention comprises an NIR spectrometer receiving at least one 
signal from a feed or product line in said process and being coupled to a computer to 
eftect continuous measurement of the spectra of the feed and/or product and provide 

25 feed back or feed forward control of the process The present invention also provides 
a computer implemented method for a system including a spectrometer linked to a 
process line containing a material X, a computer linked to the spectrometer, and a 
controller linked to the computer and the pmcess line, the computer including 
databanks having stored therein signals indicative of absorptions of standard materials 

30 (or mathematical functions thereof) and corresponding properties of said materials or 
products of said process for which X is a feed, or yield of said process, the method 
comprises steps of 

measuring absorption at more than one wavelength in the region 600-2600nm 
^ at the process line and producing absorption signals (or mathematical functions eg 
35 derivatives iherjeoO by the spectrometer in accordance therewith. 
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accessing the databanks of the computer in accordance with the absorption 
signals (or functions thereof); 

comparing, by the computer, the absorption signals (or flinciions thereoO to 
the signals (or functions thereof) of the standard materials stored in the databanks; 
5 choosing at least one standard based on the comparing, said standard having 

the smallest average value of the absolute difference at each wavelength i between the 
signal for the absorption (or function thereof) for the material and the signal (or 
function thereoO for the standards, with averaging of said properties or yields when 
more than one standard is chosen; and 

1 0 controlling said process in accordance with the outputted propeny/yield. 

The benefits of invention allow improvements in modelling with the following 
areas, identification and classification of novel products, simultaneous estimation of all 
of P properties on a sample without the need for generating P different models, and 
with the option of automatic upgrading of the model, the method being self learning or 

1 5 adjusting. The method of the invention overcomes the difficulties with the classical 
regressional approach, in particular avoiding all difficulties with numerical stability of 
the models, allowing easy and rapid identification and classification of a sample of a 
product analyzed by spectral recognition and then instant conclusions as to whether the 
sample is known or unknown, allowing simultaneous determination of many 

20 properties. In the case of blends the method can be used to determine whether the 

property of the blend is simply additive or synergetic in relation to a blend components; 
the latter is particularly useful for ditVerent blend indices and the indices considered 

The method also allows an extension of the field of application of the method 
without the need to rewrite the model, apart from the need to integrate the new 

25 samples which are inside or outside the previous field of validity of the method. This 
possibility of automatic learning, which is not possessed by traditional regression 
techniques, is a decisive advantage in the framework of continuous inline industrial 
control processes, because it allows the return of the industrial plant operations to the 
model in a certain and rapid manner in a minimum time and with all the properties 

30 considered in the model. In conirast clas.sical regression methods would necessitate 

the redevelopment of all the models, which is long and laborious without being able to 
guarantee the result of the new model obtained, because a new validation period is 
necessary, in addition during the redevelopment of the model any commercial use e.g. 
ii;i a refinery of the-model is very limited Furthermore, the method of invention allows 

35 equally the easy extension to a number of properties, which are simply incorporated 
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into the known bank. 

This remarkable possibility is tine not only for conventional properties such as 
physical chemical and/or rheolocical properties, but also for complex ones (such as 
octane number) Also it is possible to quantify by the process the response or 
5 susceptibility to lead of automobile fuels as well as the response to additives such as 
nitrates, of fuels used in diesel engines The methods of the invention equally allow 
application of the models from one apparatus to another and from one spectral region 
to another, where conventional regressive method cannot give satisfactory solutions. 
This apparatus portability is made possible by the fact that the differences between 

1 0 different spectra are the same in one apparatus as another, for the same type of 
spectrometer being considered (e.g. network scatter, Fourier transform, accousto 
optical system AOTS, diode array etc) This portability between spectral regions 
depends on the fact that as the spectral regions are intercorrelated. the relations 
between the spectra are maintained between one another. 

1 5 The invention is illustrated in the accompanying Figures in which 

Figure 1 represents a schematic diagram showing apparatus for use in the 
invention; 

Figure 2 represents a schematic block flow diagram for the method of the 
invention. 

20 In Figure I, an optical fibre links a spectrometer 2 and a probe 6 in or at 

process line I The spectrophotometer 2 produces absorbance signals at more than 1 
wavelength, which signals as such (or alter mathematical treatment to form e g 
derivative signals) are passed via line 4 to computer 5. where the signals as such or 
after conversion e.g. to one or more dei ivaiive signals, are used to enable the computer 

25 to access the databank 7 of standard signals eg absorptions and properties/yields 

therein. The signals are compared to those of one or more standard absorption(s) as 
described above and its/their corresponding property(ies) or yield(s). The output of 
the computer 5 is in the form of a property of the material in line 1 or yield or property 
of product of the process from that material and may be printed in hard copy. 

30 Preferably however, the output property/yield is used to control the process involved 
with the product in line I , ie for which line 1 is a feed or a product line; in this case the 
computer 5 is linked to and instructs the controller 9 which, via 10. controls that 
process eg. via valves/temperature and/or pressure controls in line I or in relation to 
, line 1. By this means the property of material in line 1 or yield or property of product 

35 of the process from that material can be optimised. 

22 
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In Figure 2, the initial operation 11 is to measure the absorption of the 
unknown, after which in the second step 12, the absorptions are compared to 
absorptions in spectra of standards, and in the third step 13. the spectra of the 
standards Sm are chosen according to criteria described above, and then in step 14, the 
5 property(ies) of the standard(s) Sin chosen is used to obtain the desired property or 
yield. If the spectrum of only 1 standard Sm is chosen, then the value of the 
unknown is the same as that of that standard Pm. If more than 1 spectrum is chosen, 
the value Px of the unknown is the average of the values Pm of the standards. If 
desired in an optional step 16. the value P^ is compared to the desired value for the 

1 0 unknown and in step 16 the process involving the unknown is adjusted to make the 
value Px the same as the desired value 

The invention is illustrated in the following E.xamples in which the Minimal 
Index is calculated according to the Minimal Inde.x Procedure described above. 
Mathematically the steps concerned are as follows 

1 5 For each couple of standard samples i. j, the Proximity Index ijj is determined 

from the NIR spectra by use of equation I, 2. or 3 and the properties are measured. 
For each Proximity Index is calculated the absolute ditTerence EPjj between the 
properties of the samples. The Minimal Index for property P is obtained from the 
average (EMpL) of EPjj for different values of L when L >ij. Thus the EMp(L) = 1/K 

20 Z, Zj EPjj for each of K samples for which ijj <L 

EMp(L)+ta(M) is plotted against the proximity index and in addition there is 
plotted the reproducibility of the standard method at a given level of confidence, as 
defmed in the Minimal Index Procedure above The intercept of the curve from EMpL 
and the reproducibility give the upper limit i e the Minimal Index 

25 For the Examples the data is expressed in Tables in a fonn as shown below in 

For each Proximity Index is calculated the absolute difference EPjj between the which 
the data is as follows 

30 
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The wavelengths chosen aie sliown in columns 1 and 2 
Column ?i gives the weight loading associated with each wavelength for the 
5 proximity index for the standards: I denotes no loading. 

Column 4 shows for the unknown sample the absorption at the various 

wavelengths and at the bottom the properties of that sample determined by standard 

methods. 

Column 5 shows for the unknown sample the estimated values of the properties 
1 0 and the absorptions using the method of the invention based on the properties and 
absorptions of the chosen standards 

Columns 6, 7 etc show the values of the absorptions and properties for the 
standards chosen from the bank Line 2 gives the value of the proximity index between 
the unknown sample and each of the chosen standards. 

15 
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Example 1 

Determination of Octane Number and other Properties of a Motor Fuel 

The NIR spectra between 4800 and 4000cm-' of a superfuel ID and a number 
of standard superftiels of known properties were measured. The base line was taken at 
5 4780cm- ' though similar results would be obtained with baseline drawn between 2 or 
more points. The absorbances were normalized . 

By the Minimal index Procedure described above, with use of equation 2 and 
non weighting of the absorbences the Minimal Index (Ml) was calculated to be 1 x 10" 
4. Following reference to the bank of data on superfuels and use of Procedure 1. 3 

1 0 standard samples were found with a proximity index with respect to the superfuel of 
less than Ml. The properties of these standards are shown in Table I From the 
properties of the standard samples, octane numbers (RON and MON), vapour pressure 
(hpa) volatility, percentage distilled at 70°C and at I00°C, gum content (in mg/ml), and 
content of sulphur, benzene (vol%) and MTBE were calculated for the superfuel by 

1 5 taking the arithmetic mean of the values for the 3 chosen standards. The estimated 
results are compared with the measured results. 

All the properties were obtHined from the single NIR measurement on the 
unknown superfuel and without any regression calculations, and with an accuracy in 
agreement with the reproducibilities ol'ihe reference methods Other properties can be 

20 determined in a similar way 
Example 2 

(^) Production of an Unleaded Mixed Fuel from 6 Comnonents 

A target SUPER98 superfuel olihe properties given in column 3 of Table 2al. 
was to be obtained by mixing the remain.s of a tank of finished gasoline with 5 
25 components, butane, hydrogenated sieamcracked gasoline HEN. isomerate ISOM. 
reformate (REF) and MTBE NIR absorptions at 4800-4000cm-' measured with a 
Fourier Transform spectrometer were measured, with a base line taken at 4780cm" ' 
and absorbances normalized. Results are in Table 2a 1 

Mathematic calculations were done with a computer to mix the spectra and 
30 properties of the 6 components to reproduce a finished product 

5% MTBE (on target fuel) (i.e. 4.76% in the final mixture) was "added" 
mathematically to a spectrum of the target fuel to give a mixture whose NIR spectrum 
was noted The Minimal Index was I x IQ-"* determined as described above from the 
finished gasoline 3 standards 2A. 2B and 2C were found with proximity indices with 
35 respect to the mrlicture. without weighting, and hence by averaging the properties of the 

25 
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Standards the propenies ot the mixtuie were obtained. Table 2a 2 shows the spectrum 
of the mixture, the 3 standards and the estnnation for the mixture as well as the 
properties of the standards and the estimated figures). The process was repeated with 
addition of each of the other 4 components to the spectrum of gasoline target. 
5 On the basis of the figures obtained, the blending index for each property was 

found according to the linear formula 

IP (mix) = [{ l+a)xP{mix)-P(ref)3/a 
where IP (mix) is the blending index for the ingredient in the mixture in relation to 
property P 

10 a is the percentage of ingredient in the mixture 

P (mix) is the propeny of the mixture (ingredient + gasoline) added) estimated 
by the process 

P (ret) is the propeny of the reference target gasoline. 
The blending index for addition of MTBE is shown in Table 2a3. 
5 In order to obey the linearity law here, it is necessary to limit the additions to 

not more than a quarter of the minimum to maximum range of the constituent studied 
in the industrial mixtures. However for concentration less than 20% such as for these 
oxygenated compounds, addition of 5% is acceptable. 

The process with MTBE added to the gasoline was repeated with the other 4 
20 components (and on the basis of linearity in the blending as with MTBE) to obtain 

blending indices for them as well (see Table 2a3) Then with the blending indices for 
each property for each ingredient, one can calculate the relative volume fractions 
needed to give the desired propenies for the Superfuel 9iJ and hence the blending 
order. The 6 components were then mixed in the desired proportions and then 
25 properties of the mixtui e tested and compared to those estimated by the method of the 
invention from the components present (see Results in Table 2a 4) In the estimation 
of the products and the comparison with the bank of standards, the Minimum Index 
was I x 10-"*. 3 standards 2D, 2E, 2F were found with suitable proximity indices from 
which the propenies of the superfuel were estimated by averaging as described in 
30 Procedure 1 There was good agreement between the properties obtained via the 

blending order, those measured on the fuel made and those estimated by the method of 
the invention The differences are veiy small and in the area of reproducibility of the 
standard methods. 

35 
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Example 2b 

Production of a 5 componeni leaded Supert'iiel mixture 

A target superfuel of the SUPER 97 type with 0. 15g/I of lead tetraethyl and 
having as specification an RON of 97, had an NIR spectrum as in Col 3 of Table 2bl 
5 below and other properties as given in col 3 of Table 2b I below. There were available 

4 components (HEN, ISOM, REF and an FCC cat cracker gasoline) and the remains of 
a tank of finished refined gasoline for making the target fuel. The NIR spectra of these 

5 components were measured as in Ex 2 The results are in Table 2b 1 

As in Ex.2a, mathematical calculations were done with a computer to obtain 
1 0 the spectra and properties of 5 components to reproduce a finished product. Proximity 
indices with respect to standard samples were calculated based on normalized 
absorbencies which were not weighted The method of the invention was used to find 
appropriate standards . using the procedm e of artificial mixtures as described in 
Procedure 3 and equation 8 above in wliich v was 1 and with a Min Index of 2 x 10"^. 
1 5 the latter having been calculated for standard fuel mixtures as described above 

Table 2b2 describes the results of addition of 5% of the FCC gasoline to a 
reference Super 97 gasoline target as well as the 3 standards 2G, 2H, 2J found by the 
method of this invention, from whicli the estimated properties were found. The same 
procedure was performed with the other components 
20 The blending indices were found in the same way as for Ex 2a, with the results 

for FCC gasoline in Table 2b3 and for the otlier components in the same way The 
spectral blending index (for the linear area) is obtained for each property as shown in 
Table 2b3 A blending order was also calculated, as in Ex 2a. the results being in 
Table 2b3 

25 The process of Ex 2a was repeated but with the above components and a 

Minimal Index of 2 x 10"^ The results are in Table 2b4 

3 standards 2K, 2L. 2M were found with appropriate proximity indices, which allowed 
the properties of the product to be estimated by averaging Again good agreements is 
seen between the properties estimated from the blending order and those measured on 
30 the product made, and also between the same propenies measured and those measured 
by the process The difterences seen are very small and in the area of reproducibility of 
the standard methods. Other properties can he obtained in a similar way 



35 
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Example 3 

Determination of cetane index and other properties of a uas oil 

The properties of an unknown gas oil 3 A were desired The method of this 
invention was applied with respect to a bank of known standard gas oils with known 
5 NIR spectra. The NIR spectra were obtained by F T spectrometer in the 4800 - 4000 
cm*' region [with 4780 cm"' baseline and were normalized] The proximity indices 
were calculated on the basis of Equation 2, and the Minimal Index was 2 .5 x 10"^ 
(estimated from standard gas oil data as described above). The bank of standards was 
sufficiently dense for there to be found 2 standards 3B and 3C inside the sphere with 

1 0 proximity index less than 2.5 x I O"*^ Table 3 I gives the details of the spectra and 
properties of the unknown oil A. and the standards and the estimated spectrum and 
properties, obtained by averaging. All the properties were obtained with an accuracy 
in agreement with the limits of reproducibility of the reference methods. Other 
properties can be obtained in a similar way 

15 Example 4 

On line prediction, based on NIR spectra on a mixture of crude oils fed to an 
atmosphere distillation unit, of yields iiiul properties of the diftereni distillation cuts 
such as gasoline (3S-Q5°C) benzme (^5-149=0 naphtha (149-1 75°C. jet tliel (175- 
2j2°C) light gas oil (232-242°C) heavy gas oil (342-369°C) and atmospheric residue 

20 (bp).369°C). 

An atmospheric distillation unit in a refmeiy was fed with a charge 4C which 
was a mixture in wt% of the following crudes. RUMASHKINO 8 1%, Iranian Heavy 
18%. Iranian light 1% 

Yields of various distillation cuts were desired, the boiling ranges being given 

25 above, as well as key properties of each cut as described in Table 4 .1. NIR spectra 
were measured as in Ex 1 on the crude oil Min. Index was determined from NIR 
spectra on standard crude oil (as described above) and was 2.6x 10"^ The method of 
the invention was applied using Procedure 3 and equation 8, in which v was 1. to the 
bank which was sufficiently dense for 2 standards 4A and 4B to be found with small 

30 enough proximity indices. These standards contained (wt%) (for 4A) Romashkino 
52% Iranian Heavy 29%. Arabian Heavy 1 1%. Kuwait 4%. Arabian light 2% and 
Iranian light 2%) and (for 4B) Iranian Heavy 7S%. Romashkino 21% and Arabian 
Heavy 1%. The data in Table 4 I shows the observed properties as well as the yields 
of the cuts and their properties The results obtained by this procedure were extremely 

35 satisfactory, the ditVerences observed being in accordance with standard methods of 
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measurement. Other propeaies can be obtained in a similar way 

The yields and propenies of" the distillation cuts remarkably were obtained 
directly on the basis of the NIR spectra of the feed and in line without regressional type 
calculations. 
5 Example 5 

Determination in line of the properties of a mixture of crude oils , 

Other propenies of the charye mixture of crude oils of Ex 4 were sought, based 

on the NIR spectra determined as in Ex 4 The method of the invention was applied as 

in Ex I with the Minimal Index in alt cases being 2 6 x lO"*^. Two standard crude oils 
10 4 A and 4B were found in the bank by ustnii Equation 2 The results are shown in 

Table 5. 1 . Other properties can be obtained in a similar way 

Here too the method demonstrates its capacity to predict all types of properties 

without any regression type of calculation requiring fastidious calculations. The results 

generally, as in the other Example, were in accordance with the results obtained by the 
1 5 reference methods, the deviations being found in the limits of reproducibility of the 

same methods. 

Example 6 

Determination of the Propenies of a feed to a reformer 

A feed 6D to a reformer unit was analysed by the method of the invention as 

20 described in Example 1 with the NIR spectra recorded at 2000-2500 nm. the 

absorbancies normalised and not weighted. The NIR spectaim was compared by the 
method of Procedure 3 and equation 8 (wherein v is I) with a Minimum Indexof 
2 X 10"^, which had been previously calculated as described abvoe from NIR spectra 
on standard reformer feeds Three standards OA. GB and 6C from the reference feed 

25 bank were found with small enough proximity indices; details of the spectra of the feed 
and the standards are given in Table 0. 1 . together with 5 propenies estimated for the 
feed by averaging the corresponding values of those standards The actual properties 
of the feed were measured for comparison; the measurements were by traditional 
methods (gas chromatography and density), the former necessitating laboratory 

30 determination for several hours, compared to the present NIR process which gave the 
same results in a few minutes and on line (real time in the unit) and with better 
reproducibility. 

The process allows the obtaining of a result with remarkable economy while 
(avoiding having to produce 5 regressive models The differences between the 5 
35 properties as estimated and as measured experimentally are in agreement with the 
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reproducibility of the known reference methods, namely 1.5% for jias chromatography 
for chemical compositions and 2%o for density The method can be equally applied 
for other properties such as ASTM distillation temperature curve for the feed. 
Example 7 

5 Determination of the properties of a feed to an FCC unit, as well as the yield and 
properties of the products obtained 

The NIR spectrum of the above feed 7D was measured at 4800-4000cm*^, 
with base line at 4780 cm'', normalisation of the spectrum and no weighting. The 
procedure 3 was used with equation 8, with v = I, and the Min. Index of 2.5 x 10"^ 
1 0 the latter having been previously calculated as described above from NIR spectra on 
standard FCC feeds of known properties 

The properties of the feed charge 7D sought were listed in Table 7 1 and 
included factors characterising the charge to the FCC unit, such as KUOP. crackability 
and cokability. The KUOP or Watson factor is detined as 
15 KUOP = Ve/density 60/00 

where 9 is boiling point on a Rankin scale (Absolute Fahrenheit scale) and density 
60/60 is the density of the feed at 60°F compared to that of water at 60°F 

The cracking unit operated under the following conditions riser inlet 
temperature 250°C, riser outlet temperature 525°C. MHSV (Mass Hourly Space 
20 Velocity) 78 kg/h per kg. C/0 ratio 6 6. activity of catalyst 65 (in Microactivity Test), 
The cracking gave a gasoline cut defined by ASTM distillation with initial point 
of 38°C and 90% distilled at I90°C and a residue defined by ASTM distillation with 
10% distilling at 385-C. 

By application of Procedure 3 to the bank of samples of FCC feeds 2 standards 
25 were found namely 7A. 7B and the properties and yields estimated as shown in Table 
7. 1 The results were all in line with the accuiacy based on the reference methods, as 
well as in line with the properties and yields actually ineansured Other propenies of 
the charge or products can be estimated in a similar way 
Example 8 

30 On line determination of propenies of the feed to a gasoline Hvdrouenation unit 

The gasolines obtained from steam cracking units have the inconvenience of 
containing non negligible amounts of unsaturated dienic compounds, which have the 
effect of inducing and encouraging formation of gums which are undesirable in motor 
fliel. These gasolines are therefore selectively hydrosenated to eliminate the dienes 

35 without at the same time hydrogenating other unsaturated compounds present in the 
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gasoline such as monoolefms and aromatic^ The control over these dienes is therefore 
essential not only for the final quality of the fuel (principly RON and MON) but also 
for the hydrogen consumption of the hydi ogenaiion unit. 

Units for Hydrogenating gasolines from sieamcrackers are generally coupled to 
5 a downstream distillation unit to separate a gasoline from a light cut (95% distillation 
by about 75°C) and one from a heavy cut (initial point about 95°C), before extraction 
of the benzene in the core cut and. recycle of the extraction residue from that cut called 
raffinate. 

It was desired to determine by the process of the invention the properties of the 

1 0 gasoline from the steam cracker, which was a feed to a gasoline hydrogenation unit. 
NIR spectra were obtained on the feed on line at 1000-1600 nm using a scatter 
dispersion spectrometer The absorbences were normallised. but the data was not 
weighted for use in Equation 8, in which v was I and 'min was 2.5 x 10-5 (the latter 
having been determined from NJR spectra on similar feeds of known properties). 5 

1 5 standards 8A-8E were found in the search using Procedure 3. and the properties of the 
feed calculated therefrom by averaging were all in agreement with the measured 
properties of the feed The results are shown in Table Ji. I 

In addition the chemical composiuon of the feed was obtained with great 
particularity allowing a distinction to he made for example between cyclic and non 

20 cyclic olefins as well as benzene and mono and di substituted aromatics. Equally by 
the process potential yields were obtained of the distillation cuts after the selective 
hydrogenation of the gasoline All the properties were obtained with great accuracy 
within the limits of the experimental reproducibility for that kind of property 

Other properties can be determined such as Octane Indices for the different 

25 cuts or temperatures of ASTM distillation cuives for the gasoline 
Example 9 

Method for u«; e when the density of standards in the bank i.s insutVicient 

The MON level for a reformaie 9 A was sought The NIR spectrum was 
measured at 4800-4000 cm" ' with a base line at 4780 cm" ' ; the spectra were 
30 normalised With reference to NIR spectra on reformates of known properties the 
Minimal Index was found by calculation as described above to be 2x10'^ The 
proximity indices of reformate OA and known standards were determined by Procedure 
1 The results were as given in Table I standards 9B-F were found from the 
' reformate bank with proximity indices low in relation to the reformate 9 A, but 
35 insufficiently lofv to be less than Minimal Index, as the density of the bank was too 
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small. It was thus not possible to calculate the properties with the accuracy desired. 
Procedure 1 using Equation I was replaced by Procedure 2 using Equations 4-7, with 
in Equation 4 values of Cj between -0 3 and +! .'>, in order to increase the density of 
"standards" in the bank by providing new synthetic mixtures 
5 Tables 9.2 and 9.3 show the results obtained, showing in Column 3 the 

absorbancies and properties for the "standards" (MCI, MC2) obtained by this 
densification, and with small enough proximity indices Col 4 and subsequent columns 
give the absorbances propenies of the standards 9B, 9D and 9G in the reformate bank 
used to generate the new "standards" Line 2 in these Tables show for each standard 

1 0 the fraction retained in the mixture to generate the new "standards" This fraction can 
be negative, but comprises between -0 3 (or -0 4) and + 1 .3 (in Eq.4). 

Using the data on MC I and MC2 as "standards", the properties of the 
reformate 9A were calculated by averaging (as shown in Table 9 4) The calculated 
MON of reformate 9A accords well with the experimentally measured figure, and is 

1 5 inside the limits of reproducibility of the standard method. The process can be used in 
a similar way for other properties. 

The method of the invention equally allows immediate automatic upgrading of 
the bank by automatic integration of the new samples. The process with the proximity 
indices allows consideration as a standard of all the novel "standards" introduced into 

20 the bank. This property is remarkable becau.se it allows very rapid determination of 
properties in the case of novel products nui recognised in the bank and then the 
gaining of precious time in the adjusinient of operating conditions for the 
manufacturing unit. 

Table 9.5 shows that a novel "sample" measured immediately after 

25 incorporation of the above unrecognised sample 9A in the databank, now used as a 
standard, is recognised and is perfectly calculated for the totality of its properties and 
without any modification nor intervention on the used models. It is important to note 
the superiority of the procedure over classical regressional models The latter are 
incapable of predicting propenies of samples not included within their application 

30 range or predict them with a non acceptable error, and therefore would need to be 
reactivated by the necessity to remake the model (one for each property) and this 
without guarantee of success, and with the commercial plant functioning blind during 
the recalibration period 
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TABLE 1.1 

Determination of Octane Indices anil other Properties in an automobile fuel 



Proximity lnd«x 



0 0000?7626 0.000067452 0.000072577 0.0OOO96B07 



4720 
4670 
4640 
4615 
4585 
4485 
4460 
43B5 
4332 
4305 
4260 



Benzene %Voi 
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TABLE 2:».l 

NIR Specirn of lJijJe;i(ie(l mixed fuel ;nid base fuel and additives 



0.0013833 



46tS 
4SBS 
4486 
4460 

438S 
4332 
4306 
4260 
4210 
4170 



9.0033767 



V»pow Pnttun 



4SB5 
4465 
4460 
4385 
4332 
4305 



TABLE 2a. 2 
Effect of ndditioii of MTBE on the Super Fuel 



0.00001 aoc9 



0.000060618 C 
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TABLE 2:1.3 
Blending Indices and Blending Order 



TABLE 2a.4 

Comparison of the result obtattied via the blending order and those of the 
i)ruiiiict obtained 



4616 

4566 
4485 
4460 
4385 
4332 
4305 
4260 
4210 



0 000027526 



0.000067452 



0.000072577 



0.000096607 



(blendino orderl 
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TABLE 2b.l 

NMR spectrum of Target Super Fuel, and base stocks available 
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0.OO92B28 
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0.0020427 
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0.0068411 


0.012233 


4460 


2242 


0.014409 


0.014075 


0016828 


0018119 


0.0102 


0.017615 




2281 


0 093329 


0 093268 


0095498 


0 078888 


0.095201 


O.OBS6S5 




2308 


0 14045 


0 1408 


0 1543 


0 092697 


0.14681 


0.1 1852 


430S 


2323 


0 1209G 


0 1 2085 


0 13213 


0 091686 


0.1 1947 


0 1 1078 


4260 


2347 


0.11073 


0 1 ior.2 


0 11559 


0 086428 


0.11763 


0.093662 


4?t0 


2375 


0.068913 


0 0G9012 


0 073406 


0 052861 


0.072136 


0.O6O461 


4 


70 


2398 


0 069683 


0 069948 


0073798 


0.049022 


0.076191 


0.O5776 




35 


24 le 


0 071 ;!27 


0 07W3G 


0 067143 


0.045477 


0.086471 


0.05269 




05 


2436 


0.07O00.H 


(1 070819 


0.0C5OO3 


0 050963 


0.082203 


0.056176 




2463 


0 10201 


0 10186 


0 089357 


0 13909 


0.099705 


0.1243 


4040 


2476 


0 084489 


0 08458 


0 0697 


0 1 286 


0 079381 




















RON cl«» 




97 9 


















86 2 












V 


•oour Pretiuce 




596 












Volatility 




900 4 












%Oitt lOO'c 




62.54 












%Oi«t 70*c 




42.24 
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TABLE 2b.2 

EiTect orndditioti 015% FCC gnsoline on Super 97 Product 







12 








Li 








Acuifll 


E»tim«ted 










Proximity Index 








^ "'^ °^ 


5 77E 05 

: : 


6 S3P on 

: LSL 


7 — PC nc 

iS-L 


WavalarKitfi 
















>. (nm» 















4720 


2119 




0 001495 


0 001416 


' 0 001411 






4670 


2141 




0 009125 


0.008837 


0 0O86S8 


0 008768 


0 009085 


4640 


2155 




0 009144 


0 008946 


0 009016 


0.OO8955 


0 008874 


461S 


2167 





0 0)3046 


0 012691 


0 012631_ 


0.013035 


0 01 2407 








0 009728 


0 009605 


0 009454 


0 009682 


0.009679 


4485 


2230 


1 L 


0 010956 


0 010851 


0.010639 


0.010909 


0 01 1005 


4460 


2242 


, L 


0 014524 


0 014847 


° 


0.014917 


0.0147S4 


4365 


22B1 


1 




0 09481 6 


0 094259 


0.096722 


0 093466 


4332 








0 14368 








4306 


2323 




0 1 2149 


0 1 2506 


0 12261 


0.12658 


0 12598 


4260 


2347 




0 11 096 


0 11021 




0 11141 


0.10805 


4210 


2375 




0 069127 


0 06648 


0 066676 


0 065499 


O 067266 


4170 


2398 




0 069879 




0 067881 


0 069493 


0 069464 


4135 


2418 




0 071032 


0 070704 


0 069102 


0 072032 


0 070979 


4105 


2436 




0 0G97G6 


0 06901 2 


0 07161 1 


0 069153 


0 068074 


4060 


2463 




0 1014 


0 1015 


0 10O61 


0 101 


0 10289 


4040 


2475 




0 083785 


0 081801 


0 082343 


0 082854 


0.080207 






Ratar«itce 

SUPER97 












RON ctoar 




97 a 




97 G 


97 7 


97 5 


97 6 


MON claar 




eo.j 




b6 9 


86 1 


86 7 


85.9 


Vapouf Pratiure 




6»G 




5BG 1 


590 0 


584 2 


584 2 


VotatiUtv 




906.4 




892.2 


901 0 


887 5 


888.0 


%Diit 100'e 




62.54 




62 3 


62.8 


62.6 


61.4 


%Oiit 70»e 




42.24 




41 7 


42 5 


41 7 


40.9 



In this Table 3.41 E-05 means 3.41x10-^ 



37 



wo 96/1 1399 



PCT/GB9S/02353 



TABLE 2b.3 
Bleiuliiifi! Indices and Blending order 



Volum« Friction 



87.9 
235. 0 
204.7 



TABLE 2b.4 

Comparison between the resnlts from the blending order and the product 
obtained 







We.jiht _ 




2K 


2L 


2M 








Meatiirert 


Estimaiea 








ProximitY Indax 








8 41E 05 


5 35€-05 


0.000159 


0.000199 


















X (em- 1 1 
















47 20 


2119 




0 001742 


0 001809 


0.001776 


0.001851 


0.001802 


4670 


2141 




0 009508 


0 009166 


0 009131 


0 009262 


0 009104 


4640 


2155 




0 008G98 


0 00926 


0 0O9326 


0.009695 


0 008759 


4616 
4585 


2167 




0.012758 


0 013206 


0 012772 


0.013574 


0.01327 


4485 


2181 
2230 




0 009725 
0 010469 


0 009649 
0 010438 


0 00924 2 
0 010486 


0 009921 
0 010336 


0 009484 
0 010493 


4460 


2242 




0014142 


0 014252 


0 014878 


0 013896 


0 013982 


4385 


2281 




0 090899 


0 093317 


0 092897 


0.09338 


0093676 


4332 


2 308 




0 1 3680 


0 13974 


0 1 3652 


0 1 3569 


0 14683 


430S 


2323 




0 1 159G 


0 12103 


0 1 195 


0 1 248 


0 11879 


4260 


2347 




0 1 1499 


0 1 12H 


0 11271 


0 11216 


0 11352 


4210 


237G 




0 071524 


0 0G8713 


0 071255 


0 068318 


0 066566 


4170 


2398 




0 07OCC2 


0 070304 


0 069855 


0 069368 


0 071669 


4135 






0 072077 


0 071069 


0 070001 


0 07464 


0 068566 


4105 


24 36 




0 0C944B 


0 07175G 


0 073371 


0 071076 


0 070822 


4060 


24G3 




0 10444 


0 10095 


0 10163 


0 098515 


0 10269 


4040 


2475 




0 OBGI 16 


0 08264 


0 084647 


0 083297 


0 079975 






















st.lndards 










RON Clear 


97 7 


97 5 


97 7 


97 8 


97 5 


97.7 


MON clear 


UG 3 


8C4 


8G 2 


86 


86 1 


865 


Vapour Preature 


S89 9 


595 


598 0 


596 


600 


598 


Volatilrty 


952 7 


949 


956 7 


955 


960 


955 


%0..t 100'c 


65 5 


62 


63 0 


66 


63 


60 


%Oiat 70«e 


45 4 


47 


44 2 




466 
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TABLE 3.1 

Deteriniiintion ofcetane index ami other properties of a gas oil 









Gat 


Oil A 


3B 


3C 




Proirimitv Index 




Measured 


Eitimaied 

1.71E-C6 


1.39E-06 


2.23E-06 












• 




X (cm- 11 


X (nml 












4720 


2118.6 


1 


0 000120383 


0 000126616 


0.000139825 


0.0001 1341 1 


4672 


2140.4 


1 


0 001962853 


0 002013913 


0.002015676 


0.0O2O1 1949 


4640 


2166.2 


1 


0 0034 34 74 7 


0 003416109 


0.003438675 




4616 


2166.4 


] 


0 004544314 


0.OO4490799 




0.0O45OS037 


4584 


2161.6 


1 


0 004729896 


0 0046754 


0.00463465 




44S4 


2230.2 


1 


0 0071 19883 


0.0O6932337 


0.006908771 


0.0O69559O3 


4460 


2242.2 


1 


0 010349409 


0 010133388 


0 010064653 


0.010202122 


4364 


2281 




0 0 74C06084 


0 074925207 


0 0749301 1 7 


0074930298 


4332 


2308 4 




0 16BG7785? 


0 157745031 


0 15799051 


0.1S7499551 




2323 4 




0 101624835 


0 102266697 


0.102217602 


0 10231S793 


4260 


2347 4 




0 1 J1871507 


0 131331453 


0.131380548 


0 131282357 


420B 


2376.4 




0.088627865 


0 088637684 


0 08855913 


0.088716237 


4168 


2399.2 




0 092B99205 


0 093012126 


0.093105408 


0.092918844 


4132 


2420.1 




0 084003812 


0 084 7691 1 


0.084827844 


0.084690376 


4104 


2436.6 




0.081106377 


0.081 2978S1 


0.081361676 


0.081234027 


4060 


2463. 1 




0 08642837 


0.086487285 


0.086487285 


0.086487285 


4O40 


2475.2 




0 007 192608 


0.O6740372 


0.067418449 


0.067388991 


Cetsne Index 




52 


52.5 


53.3 


51.7 


Cai*ne Number 




55.3 


52.75 


52.1 


53.4 


Oeneitv 15-C 




0 8434 


0.84085 


0.8385 


0 8432 


Fta.h Point 




62 


57.6 


60 


65 


% Sulphur 




0.29 


0.25 


0.23 


0.27 


Cloud Point 




5.1 


5.6 


5 


6 


Filtiabilitv 






0.5 




O 


Vi«cotitv 40*C 




3.1 


3.7 


3.7 


3.7 
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TABLE 4.1 

Determination ot'yields :ind properties ui cuts troni distillation or mixture or 



crude feed oils 







We,,ir,t 


Ctiarte 4C 












Meaturad 


Estimated 


4A 


48 




ProximHY Inrtex 






9.98E07 


1.21E-06 


1 33E-06 
















k Icm-ll 


>. (nml 












4672 


2140.4 


1 


0.001777942 


0.001748627 


0.001771733 


0.00172562 


4640 


21S5.2 


1 


0.003139917 


0.003211964 


0.003256211 


0.003167717 


4616 


2166.4 


1 


0.00377911 


0.003827795 


0.00383S639 


0.0O38199S2 


4584 


2181.6 


1 


0 003794844 


0.003797791 


0.003829737 


0.003765845 


4484 


2230,2 


1 


0.0O6O94969 


0.00614464 


0.006272386 


0.006016694 


4460 


2242.2 


1 


0.009268476 


0.009156818 


0.009276767 


0.009034879 


43B4 


2281 


1 


0 078089814 


0 077898738 


0.077667019 


0.078130467 


4332 


2308 4 


1 


0 15773336 


0 157794497 


0.157793411 


0.157795S84 


4304 


2323 4 


1 


0 104631 107 


0 1045241 


0 104179066 


0.104869135 


4 260 


2347 4 


1 


0 130690546 


0 130249322 


0 130445176 


0.130053468 


4208 


23764 


1 


0087815393 


0.08775 1064 


0 087838988 


0 08766312 


4172 


2396 S 


1 


0 091208037 


0 090879399 


0 090878774 


0 090880025 


4132 


24 20 1 




0 084048925 


0 084706329 


0 08465164 


0 084 761019 




243G.C 




00824855 


0 082364989 


0.082389016 


0 082340962 


4060 


2463 1 




0 087068028 


0 087578898 


0 087475 


0.087682795 


4040 


2475 2 




0.007784043 


0 068366138 


0 068439449 


0 068292827 


Oantitv »5*C 




0.8663 


0,86655 


0.8646 


0 8665 


% Gacokn« 








7.4 


7.4 






7.0 


7.2 


7.3 


7.2 


% Napht* 




4.3 


4.5 


4.5 


4.5 






8.5 


8.6 


8.6 




% liOM gas oil LGO 




18.9 


18.8 


19.2 


18.5 


% Heavy n«, oh 




4 5 


4.5 


4.6 




% Retidu« RAT 




49 


49.2 


48 6 


49.8 


% Parallinei Naphta 




52.2 


52.1 


53.1 


51.2 


FUih Hoint Petrol 




59.2 


59 5 


59.8 


59.3 


CIDUO point LOO 




•8.1 


•8.1 


•8.8 


•7,5 


% Sulphur RAT 




2.8 


2.8 


2.8 


2.9 


Viscontv 100'C HAT 




53 16 


52.72 


48 53 


56,91 
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TABLE 5.1 

On line determination of properties ot'n mixture of crude oils 











"Eiti'mated 








PfO«imHy Index 










■■ : — 






























4672 


2140.4 




0.001777942 


0 001 747949 


0 001 771 05B 




4640 


2155,2 




0.003139917 




0 003254971 


o on^ K Ktt 


4616 


2166 4 




0.00377911 






0,00361 B445 


45B4 


2181. S 




0.003794844 


0.0O3796319 


0.003828279 


0,00376436 


4484 


2230.2 




0.006094959 


0.006142161 


0 006269999 


0.006014322 


4460 


2242 2 




0.00926B476 


0 009162272 


0 009273227 


0,009031317 


4384 


2281 




0 078089814 


0.077868554 


0 077637461 


0.078O99647 


4333 


2308 4 




0 10773336 


0.1077333C 


0 15773336 


0.1S773336 


4304 


2323.4 




0 104631 107 


0 1044836 


0 104139419 


0 104827782 


4260 


2347.4 




0 1 30690546 


0 130198858 


0.130395533 


0.130002163 


4208 


2376 4 




0 OB7815393 


0 087717055 


0 087805559 


0 08762B5S2 


4T72 


2396 9 




0 091208037 


0.090844188 


0.090844188 


0.090844188 


4132 


24 20 1 




0 064048925 


0 064G7361 


0 084619424 


0.084727595 


4104 


2436 6 




0 0824855 


0 062333077 


0.082357661 


0.082308493 


4060 






0 OB70C8028 


0.087644965 


0.087441 71 


0.087648219 


4040 


2475 2 




0 067784043 


0.06833965 


0 0CB413403 


0.068265897 


Oaniity 


1 0.8663 


0 86655 


0 8646 


0.6665 


% Sulphur 


! 1.6 


1 65 


1.6 








2.27 


2.265 


2.36 


2.17 


% Conr«d»on C.vbon 




4.8 




4 8 


6,1 






5 


4.95 




5 



TABLE 6.1 

Deterniinatiun of the Properties of a reformer feed 



0 33C44352 
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TABLE 7.1 

Determinntioii of Properties ui't'eed tu FCC reactui* and yields and properties of 

products 







Weinrn 


tsi 







7B 








M«»«ured — 


Ettimaied 








f roKimily Indan 






'^Ot-06 


'^BEOfi 


1 ,30E06 
















X. (em-1) 














4720 


2118.6 




0,00024017 


0.000283004 


0.000238346 


0.000327662 


4672 


2140.4 




0.002238801 


0.002010364 


0.001890879 


0.00212985 


4640 


2155.2 




0.004237234 


0.003903227 


0.003874 1 1 7 


0.003932336 


4612 


2168.3 




0.005237444 


0 0O4372667 


0.004866233 


0.0O50791O2 


4584 


2181.6 




0 005332797 


0.005055096 


0.005031089 


0.005079102 


4464 


2230.2 




0.007970887 


0 007756354 


0.007744 263 


0.007768446 


4460 


2242.2 




O.On 303204 


0.011210967 


0.011212199 


0.011209736 


4384 


2281 




0 072994455 


0.07273491 


0.07292398 


0.072545839 


4332 


2308.4 


1 


0.152067643 


0.152159348 


0,151945649 


0.152373O47 


4304 


2323.4 


1 


0.100517606 


0.100397569 


0.10O601923 


0.100193214 


4260 


2347 4 




0 131209247 


0.131514201 


0 131487607 


0 131540794 


4212 


2374 2 




0 091618024 


0 091623192 


0.091564633 


0.091681751 


4168 


2399.2 




0 094011773 


0 094322962 


O.094275B2 


0.094370104 




2420.1 




O.OU61«4908 


0 086G75314 


0,086678538 


0.08667209 


4104 


2436.6 




0.081457005 


0 081916022 


0.081981 133 


0.081850912 


4060 


2463.1 




0.0842G7922 


0.O84318052 


O.084444O43 


0.0S4192O6 


4040 


2475.2 




0.06911082 


0 0G91467S2 


0.069239547 


0.0690S3957 






0 926 


0.9225 


0.922 


0.923 


% Sulphur 




1 97 


1.85 


1.83 


1.87 






«35 




78 2 


88.2 


Vi<co«i(y 100'C 




0 e 




8 7 


9.5 


Temp. 60% disiillett 






464 






KUOP 




1 1 8 


1 1.85 


118 


11.9 


Mol Weiohi 




4G0.C 


449 95 


434,5 


465 4 


% Aromot.c C.iibo.i 




<;i.8 


21.2 


21,6 


20.9 


CHACKABILITY 




y 4 7 


2 57 


2 55 


2,59 


COKABILITY 




1 01 


1 00 


0 99 


1,01 


GASOIL INDEX 






1 515 




1 49 


GASOLINE INDEX 




0 99 


0 985 


0 99 


0 98 


Gstoline Yield l%l 






44 5 




46 


Re«klu« Yield (%) 




12 


11,25 




9 5 


RON Clear Uasoliiie 




'J2 7 


92 4 


92.4 


92,4 
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TABLE 8.1 

On line Determination of properties oMeed to hydi ogeiiation unit for gasoline 



0.0000034 
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TABLE 8.1 (cout) 

On line Determiuntion of properties of teed to liydrogenation unit for gasoline 











" ^Ettimated 




, 




22 


Si 








Manured 








_ 








Prox (nde« 






0 OO00O4 2 


0,0000034 


0 0OO0O47 


0 OQOOl — 
: Li. 


— T 

0,00001 30 

























A. (cm-n 


X. tnml 




















8140 


1226 5 





0 00894 


0 0089636 


0 0Q9 


0,006866 


0.009078 


0,008757 


0 009014 


806S 


Tjto 





0 004327 


0 004-t40r 


0 0044fi7 


0 00427 


0.004393 


0.00421 


0.004363 


7768 


Tiii" 





0000897 


0 OQQQ-QB 


n 


0 0009^ 


0.0OO918 


0.001025 


0.000921 


8117 


T2J2 


- 


0 00534 


°0 008592 




0 0O5591 


0 005209 


0.005491 


0.005213 




TmT 




0 00869 






0,0O8798 




O.OOB701 


0.00846 


7396 


1352 




0 01 2209 


0 012095 


0 012056 


0,012119 


— °o°onB2 


0.01 2289 


0.01 2093 




1356 




0 01580G 


0 015784 


0,015674 


0015614 


0.015688 


0.016136 


0.01581 


7356 


T 359.5 




0.022613 


0 022831 


0,022633 


0.022499 


0.022916 


0.023452 


0.022666 


734B 


1361 




0.024681 


0.024929 


0024753 


0.024623 


0.025033 


0. 025638 


0.024698 


7339 


1362,6 




0,026435 


0.0267O7 


0 026562 


0.026428 


0.026818 


0.027301 


0.026438 


7321 


1366 




0.029615 


0.029773 


0 02976 


0.02954 


0.029B49 


0.030217 


0.029509 


7273 


1375 




0 038104 


0.038252 


0 038215 


0 037747 


0.038738 


0.038781 


0.037777 


7254 


1378.5 




0.O42097 


0 0423C1 


0,042333 


0 041789 


0.0429O4 


0.04293 


0.041847 


7241 






0 044261 


0 044611 


0,044521 


0.043951 


0.044955 


0.045039 


0.044087 


7231 


1383 




0.0454 1 


0 04 5560 


0 045625 


0 045014 


0.045B82 


0.045992 


0.045269 


7199 


1389 




0 04833 


0.04827G 


0,048333 


0 04766 


0.048684 


0.O4862 


0.048084 


7184 


1392 




0.049461 


0.049347 


0,049442 


0.048812 


0.04B816 


0.049597 


0.04907 


7161 


1396,5 




0,049514 


0 049304 


0,049526 


0.049149 


0.049522 


0.049216 


0.049109 


7151 


1398.6 




0 04893 


0.048706 


0,048937 


0.048708 


0.048735 


0.048454 


0.048697 


7117 


1405 




0.0471 


0,04703 


0,047083 


0.046966 


0.046646 


0.046465 


0.048003 


7106 


1407.5 




0.046793 


0 04G839 


0.046789 


0.046579 


0.046484 


0.046242 


0.O481 


7087 






0 045R55 


0 046047 


0045894 


004S586 






0.047474 
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TABLE 8.1 (cont) 

On line Determination oC prupcrties oCi'eed to hydrogeuatiou unit for gasoline 







Weipht 


feed OF 


«A aa 8C 


§2 


8e 








Ueatiiien 


Estiinaiecl 














Prox. Index 






0 000004 2 


0 0000034 


0 000004 7 


0 00001 14 


0.0000130 


0.0000222 






















>. (cm- 11 




















7070 


1414.6 




0 04 3481 


0 043G«2 


0 043G73 


0 04322 


0.043476 


0.043002 


0. 046138 


7018 


1425 




0.036493 


0 0JC481 


0 036252 


0 036396 


0 036261 


0.03S933 


O.037S64 




1430 6 




0 037461 


0.037372 


0.037024 


0.0375 


0.037441 


0.03717 


O.0377J4 


6974 


1434 




0.037514 


0.037391 


0.037162 


0.0376O5 


0.03762 


0.037269 


0.037298 


6971 


1434.5 




0037307 


0.0372G1 


0.03706 


0037486 


0.037496 


0.037136 


0.037127 


6930 


1443 




0.032307 


0 032112 


0 032218 


0.032451 


0.032217 


0.031884 


0.03179 


6649 


1460 




0.022126 


0.021897 


0 022177 


0.022295 


0.021806 


0.021642 


0.021665 


6824 


1465.6 




0 021691 


0.021377 


0 021682 


0 021745 


0 02138 


0.021154 


0.020925 


6762 


1481 




0 018013 


0 0179G 


0 018203 


0 01841 1 


0,017735 


0.017774 


0.017678 


6720 


14B8 




0.01 609B 


0 OIGOOe 


0 016257 


0 016607 


0 015759 


0 015855 


0 016811 


6G93 


1494 




0.014581 


0.014542 


0.014741 


0.01504 2 


0,01436 


0.01435 


0.014215 


6614 


1512 




0.010568 


0 01046 


0.010646 


0.011075 


0.010321 


0 010193 


0.010114 


6566 


1523 




0.007881 


0 007794 


0 007827 


0.008437 


0,007533 


0 007562 


0 007616 


6S36 


1530 




0.0071 1C 


0 007114 


0 00711 


0 007727 


0 006783 


0 006993 


0.006957 


6461 


1543 




0.00523:i 


O.0O62GG4 


0 006299 


0.005895 


0 004821 


0.005154 


0.005163 
















































•A Imeur Saturatefl 


8 73 




7 73 


7 98 


10 65 


8.86 


9.69 




6.83 


7.35 


6.51 


6 G2 


8.8 


8 26 


6.56 




6 83 


6 T.s 


C 38 


C 93 


5.56 


4.67 


4.95 


% Linear Olstiiu 


11 33 


11,92 




12 09 


10 49 


1 1.33 


13.28 


•A Cyclic OlBfiiif 


1 2 79 






13. GB 


12.94 


1 1.81 


12.44 


•A B«i<iei>e 


2:^ 4G 


21 74 


21.35 


21 38 


23.34 


21.63 


21.1 


% Tolueite 


13 72 






13.51 


13.55 


13.18 


14.18 


% Xvfene 


5 08 


G 74 


5.97 


5 57 


4.89 


1 1 36 


5.9 


% Alkvl beniene 


6 9G 


6 94 


5.98 


0 6 


5 11 


6 84 


6.16 


% Dieiies 


IG 44 


10 90 


1 7.67 


17 94 


14 83 


13.67 


15.47 


OENSITV 


0 8124 


0 8000 


0 8133 


0 8097 


0.8012 


0.803 


0.806B 


Yi«lc1 hflht CU1 


26 7 


2R 7 


24 C 


25 1 


28.6 


27.6 


27.B 


Yield heavy cut 


40 3 




41 7 


41 2 


35 7 


39.7 


38.7 




13 7 




13.5 




13 7 


12.4 


13.5 


Yield ben7ene 


:'0 3 


20 G 


JO 2 


20 3 


22.1 


20 4 


20 0 
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TABLli 9.1 
Deterininatioii ol MON of a reformnte 







W«iQht 


R«lor 
Meatiired 


nste 9A 


96 


9C 


90 


BE 













Eitimatod 
0.000058424 


0,00000398 


0.00006638 


0.0OO18149 


0.00018529 


0.00019385 


X (cm-11 


^ ^^^^ 


















4720 
4670 


2119 


- 


0,0010961 




0.0012043 


0 0012359 


0 00076051 


0.0010861 


0.00091553 


4S8S 


2155 
2181 


• 


0.018144 




0.019076 
0.025324 


0.016368 
0.021 B03 


0.014737 
0.014281 


0.01558 
0,0169 
0.023414 


0.015681 
0.015344 
0.021006 


«4eB 
4460 


2230 
2242 




0.020515 
0.012619 
0,018197 




0.020CI2 
0.011885 
0 017086 


0 017869 
0 011415 
0.016221 


0016021 
0.010242 
0015362 


0.016114 
0,01136 
0.01fi7»e 


0.016429 
0.010629 
0.016344 


438S 


2281 




0,092064 




0 089813 


0.090082 


0.09165 


0.094457 


0.088586 


4332 






0 1 2886 




0.12812 


0 13172 


0,12815 


0. 1 269 


0. 1 3463 


4305 


2323 




0 1 1882 




0 1100C 


0 11689 


0, 1 205 


0 12271 


0 11448 


4260 


2347 




0 098322 




0 0971C 


0 10023 


0 099018 


0 09946 


0.10223 


4210 


2375 




0.064577 




0.063284 


0 005124 


0 066654 


0.065915 


0.06736 


4170 


239B 




0 06140G 






0 0634 25 


0067026 


0.065132 


0 065986 


4135 


2418 




0,059290 




0.058485 


0 061 794 


0 066184 


0.063244 


006268 


4105 


2436 
2463 




0.06198 




0 061976 


0 064438 


0.066333 


0.062871 


0 066088 




2475 




0.1 1037 
0 091698 




0 11622 
0.094281 


0 0916C1 


0,11289 
0 091751 


0.10222 
0.096914 


0.1122 
0.090604 


MON 




0 


m 




H8 3 


86.2 


87 2 


89.2 


82 4 



TABLE 9.2 
"Sample" MCI obtained by densiUcation 



4585 
4485 ~ 
4460 ~ 



O 01 77G778;' I 
0 0241 1R931 
0 019O41012 



0 1300549 



0.0G0ierii90 
0 or.24r.tiGf.i> 



46 



wo 96/11399 



PCT/GB95/023S3 



TABLE y.3 
"Silinple" IMCl ubciiiivii by deiisificsuion 




1 ABLI£ V.4 

Detei iniiintioii ofMON ori crorm:Ue based on "Snmples" generated 



4720 
4670 
4640 
4615 
4585 
4485 
4460 

4332 
4305 
4260 
4210 
4170 
4)35 
4105 
4060 
4040 


21 19 

2155 
2167 
2181 
2230 
2242 
2281 
2308 
2323 
2347 
2375 
239B 
2418 
2436 
2463 
2475 




hetof 

0 0010981 
U 0) 7744 
0 018144 
0 024297 
0 020515 
0 0)2619 
0 018197 
0 092004 
1) 1 28B6 

0 098322 
CJ 004577 
() 001405 
0 05929G 
(1 0C198 
0 1 1037 
0.091698 


iial« 9A 

0 00OOO9271 

0 001 1773 
0.017321 
0 017B4G 
00239335 
00192925 
0.01 1968 
00171385 
0 0916295 
0 130795 
0 1 1671 
0 09869 7 
0 064435 
0.0617075 
0.06004 75 
0.06264 
0 110815 
0.092348 


MCI 
0 0O0O1O043 

0.0011759 
0.017366 
0 017768 
0 0241 19 
0.019041 
0 01 1986 
0 017123 
0 09183 
0. 1 3OO0 
0 1 191 2 
0 098596 
0 06452 
0,061987 
0 000186 
0 062409 
0 1 1024 
0 093416 


MC2 
0.OOOO13457 

0.0011787 
0.017276 
0.017922 
0.023748 
0.019544 
0 01 1948 
0 017154 
0 091429 

0.13153 
0.1 183 
0 098798 

0.064 35 
0.061428 
0.059909 
0 06281 1 

0 1 1 1 39 

0 091 28 


MONOO 




0 


88 


88.3 


86 4 


68.2 



wo 96/1 1399 



PCT/GB95/02353 



TABLE 9.5 

Recogiiiliuii ufiu'w sniitple by aiilu ndjiislineiil 



O 00001223S 
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Example 10 

On line Determination of the Prop erties ot'Polvbutenes in a Polvisobiitene Plant 

It is desired to control the properties ot'polybiitenes made during their 
manufacture in order to adapt immediately the operating conditions to any changes in 
any product. 

They were made by polymerisation ot isobutene with an Bronsted acid catalyst 
to form a crude product from which distillation removes gaseous hydrocarbons and 
light polymeric products, and leaves heavy polyisobutene 

During the manufacture, the absorbances of heavy polyisobutene (Ref lOA) 
were measured with an NIR spectrometer in the wavelength region 6000-8850cm- ' 
The spectrometer had been installed on line in a plant with the aid of a fast side loop 
situated in the line cari-ying the heavy polymer remaining arter the distillation. An 
analyser attached to the spectrometer sent within 2 minutes to the controller of the 
plant a measure of the propenies which had been obtained on the basis of the NIR 
spectrum measured on the above products and estimated using the method of the 
invention as follows. 

The method chosen to treat the NIR spectrum involved a discrete selection of 
wavelengths chosen on the basis of chemical and/or statistical criteria, the chosen 
wavelength being between 6070 and SSO.lcm.- ' The absorbances were normalised 
according to procedure (2) 

For a series of standard polvis^buiene products, tor which the NIR spectra 
were known, the Minimum Proximity Index was obtained by the method above to be 
5xlO-^; this Index was not weighted The proximity indices, between the absorbances 
of the standards in the bank and those of the unknown from the plant (normalised as 
above) were calculated and 5 standards lOB-lOF were found with proximity indices < 
Min. Prox Index 

By averagmg the values of each property of those 5 standards it was possible to 
calculate the corresponding properties sought for the unknown for controlling the plant 
from knowing the viscosity at !00°C. the number average molecular weight (MN). size 
of the distribution of molecular weights obtained by gel permeation chromatography 
(called LGPC) as well as the content ofbutene-l 
(BUT-1) 

Other propenies may also be determined such as density or inflammability point 
in a similar way ' 

Table 10"I shows the results, from which ii is clear that the calculated 
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properties (in col 5) were all in agreement with those measured on the unknown by 
standard methods and were within the limits ofreproducibiHty of those methods for 
use in the polyisobutene area (on the basis of a 95% probability in any measurement) 
namely 0.7% for viscosity (ASTM D445), 8 5°C for inflammability point (ASTM D93- 
5 80), 5% for number average Molecular Weight and 3% for Molecular Weight 

distribution (both by Gel Permeation Chromatography), iO% for unsaturation (NMR) 
and for buiene-1 (IR) and 1% for maleinisation index. 
Example 1 1 

Determination of Properties of Low Molecular Weight Polyisobutene 

10 The propenies of a low molecular weight polyisobutene ( i I A) were 

determined in order to control its manufacturing unit as in Example 10 The method 
adopted was as in E.xample 10 with the absorbances of the polyisobutene measured as 
before in the 6000-S850cm* ' wavelength range with the aid of the NIR spectrometer 
installed on the residue line from the distillation unit. 

1 5 The propenies for determination were the viscosity at 100°C, the number 

average molecular weight, the LGPC" (as in Ex. 10), the content of butene-1, the 
inrtammability point (IP) and the degree of unsaturation expressed in groups/litre, and 
the maleinisation index (PIBSA). The unsaturations were of the types VIN, TRII, 
TR12, TRI2cis, TRl2trans. TRITotal and TETRA as defined above The maleinisation 

20 index is particularly important for control of plants to make polyisobutenes as it is of 
great value to purchasers of low polyisobutenes. 

The NIR spectra of a series of standard polyisobutenes whose properties were 
measured by reference techniques, were determined and from the bank obtained the 
Minimal Index was determined at 9x lO"-*^. via the unweighted proximity indices. The 

25 density of standards in the Bank was sutViciently high for there to be 5 standards 1 18- 
1 1 F inside the sphere with proximity indices with respect lo the unknown less than the 
Minimal Index. By averaging the data from these standards, the properties of the 
unknown were determined The results are shown in Table 1 1 I Here also the 
properties of the polyisobutene from the plant calculated from the standards were all 

30 within the limits of reproducibility of the standard methods Thus continuously and 

with total reliability, the quality of product from the plant can thus be obtained and can 
be maintained taking account of the process dynamics. 

Other properties of the low molecular weight polyisobutene may be obtained in 
, a similar way 

35 

50 
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Example 12 

Determination of the Propenies of a Polyethylene to Control lis Manufacture 

Ethylene was polymerised with a chromium catalyst in a plant to produce 
polyethylene, whose propenies of density and tUiidity index/Melt index (measured 
5 according to the 1F2 standard method) (called grade) were to control its market 

specification. The product 12 A. in the form of a powder, required, rapid measurement 
of these properties to correct the plant operating conditions to ensure manufacture of a 
polymer of constant quality 

The NIR spectra of a series of standard polyethylenes were determined by 

1 0 means of a Fourier Transform NIR spectrometer in the 5500-8400cm- ' region, as well 
as their densities and fluidity indices to obtain a bank of standards. The absorbances of 
the spectra were normalised as described above, to ensure better known numerical 
stability for the data The Minimal Proximity Index was calculated from the 
unweighted proximity indices of the standards (according to Eci.2) and by means of the 

1 5 technique of Eq.S in which v wa.s 11. tlie pn^ximiiy index to the product I2A was 
chosen at 0.025 

The proximity indices between the unknown polyethylene 12A from the plant 
and the standards were calculated according to £c|.2. and three standards I2B, 12C 
and 12D were found with small enough proximity indices From the properties of 

20 these standards 12B-D by averaging, the properties of the unknown product 12A were 
obtained in less than I minute, allowing inmiediate reaction to all variations in the 
production operation The results are shown in Table 12 1. and are in perfect 
agreement with the propenies determmed by reference methods, and within their 
reproducibilities, namely l%(i for density and 14% for grade. 

25 The method may be applied m the same way to the determination of other 

properties, for example percentage of conionomer in an ethylene copolymerisation, the 
degree of conversion of the reaction or the content of volatiles. as well as to other 
types of polymerisations to polyethylene .such as ones with Ziegler Natta catalysts 
Example 13 

30 Control of Polyisobutene Production 

Polyisobutene production was numitorecl as described in Example I 1 over a 
period of 24 hrs with a read out thi oughoui of the calculated viscosity of the product 
heavy polymer Fig. 3 shows this reatl out from the computer analysing the NIR data 
with viscosity expressed in centistokes as ordinate and time in hours as abcissa reading 

35 from right to lefti-'together with a copy of the spectrum (with linear ordinate range 
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from 0-0 5 absoibance unit and 950 to 1650 nni abcissa measuring from left to right) 
of the product at the start of the 24hr period when its properties were as follows: 
viscosity 2538cps, Flash Pt 197°C, Mn 2267, Mol Wt Distrib. 1 93, % vinylidene 
3.364. TRITOT 60 03, TRI2C 12.63. TRI2t 36 69% butene- 1 2 65 Over this 24hr 
5 period the viscosity changed significantly The calculations to enable this to be 
followed continuously were possible by the method of this invention involving the 
remarkable property of the process of "automatic densificaiion" and auto-adjustment 
of the standards, with avoidance of the recalibrations and fastidious and uncertain 
remodellings of regressional methods. This aspect of the present method is described 
1 0 further below 

For many process operations involving on line analysis of product leaving the 
reactor, there may be only a few standards involving NIR absorbances on these 
products of known analysis or properties 

For a polyisobutene pilot plant, the absorbances of the standards for the 
1 5 products, were measured with an FTIR spectrometer between 6400-9000cm"' and the 
properties were the percentages of constituents leaving the reactor as measured by 
distillation (i.e. the heavy and light fractions); these parameters were measured because 
they represent the most interesting to follow for the progress of the reaction, though 
others could also be used as well or instead if desired The Minimal Index was 
20 calculated from the proximity indices of the standards (according to Eq 3) weighted 
here by the accuracy factor in absorpiiuns ai each wavelength chosen, the proximity 
index as with respect to the unknown wns fi.\ed by Procedure 3 at 3x10"-* Because of 
the low population of standards in the band with proximity indices with respect to an 
unknown polyisobutene I3A from the unit within this level, no standards were found 
25 and. hence no values can be directly esiiiiiaied for the properties (shown as a question 
mark in col. 5 of Table 13.1 hereafter with the rest of the absorbances of the unknown 
and those of 3 nearby standards 13B. I3C and I3D. 

The densification process of procedure 2 was thus invoked in order to increase 
the bank of known standards Tables 13 2 and 13 3 (in col 3) show the results 
30 obtained for "standards" MCI and MC!2 respectively determined by calculation and 
located with small enough proximity indices of process Eq. 3. Col 4 and later 
columns give the absorbances and properties of the existing standards of the bank used 
to generate novel "standards" namely 138, 13C! and 13E (for MCI) and 13B, I3D and 
^ 13E (for MC2). Col 2 gives for each of the novel standards the fraction from the 
35 known standards used in the mixture to generate the new ones This fraction can be 
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negative or positive but can be applied as in Eq.4. 

Based on the new "standards" MCI and MC2 (with Proximity Indices with 

respect to the unknown of 2.7x10"^ and 2.6x10-^ respectively), the properties of the 

unknown 13A were then calculated, by use of arithmetic averaging with results as 
5 shown in Table 13 4. The weight percentage of the 2 polymeric components in 

product 13 A were well compatible with the measured values from.the distillation in the 

light of the errors in the reproducibility of the latter 

The above process enables immediate auto adjustment of the bank by automatic 

incorporation of the novel "standards" as though they were originally measured ones. 
1 0 Thus property is remarkable as it enables determination of properties of products 

hitherto unknown in the bank and to gain precious time in the adjustment of operating 

conditions on the plant to take account of changes. 

This autoadjustmeni process was used in a process whose results were as given 

in Table 13.5 below, which shows thai the novel "standard" measured immediately 
1 5 after incorporation of the "standard" unrecognised in the bank is now recognised and 

can be used to calculate well all the properties and without any modification nor 

intervention from existing models. This operation is better than that of classical 

regressional models, which are incapable of predictmg the properties of samples not 

included in the trial area or of predicting them with a non acceptable error ; these 
20 regression techniques would require for use the preparation of new models for each 

property and without guarantee of success and the production unit using the new 

model would be operating blind during the recalibnition period. 

Example 14 

Determination of the Propemes of a Polyethylene Glycol 

25 A polyalkoxylenated product had been made discontinuously by polymerisation 

in the liquid phase of one or more epoxides specifically ethylene oxide with an organic 
compound possessing at least one active hydrogen atom, such as an alcohol specifically 
butanol. The values of the properties of the product during the process had hitherto 
been regularly evaluated in the laboratoi-y by standard methods during the production 

30 in order to determine the necessary' amounts of epoxide However, the times to obtain 
the analytical results were generally prohibitive in terms of the amounts of epoxide 
consumed and quality of product obtained during non optimum operations 

The method of the present invention was applied to this process A band of 
standards for polyethylenoxylated butanols covering the field between low ones (Mol. 

35 Wt about 200) up- to high ones (Mol. Wt of the order of 9000) was generated 
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incorporating the properties of a number of these compounds as well as their spectra, 
determined in the 4000-8400cnr ' wavelength region with an FT NIR spectrometer. 
The properties considered were the hydroxyl index (fundamental for the 
conduct/progress of the reaction) as well as the viscosity measured at 100°C and the 
5 molecular weight. 

The bank of data was applied by the method of the invention to the production 
of a polyethylenoxylated butanol of Molecular Weight of about 8000 (PEG 8000). 
From the standards with spectra normalised per Eq 3 the Minimal Index was calculated 
at 1.2x10"^. The Proximity Indices between the unknown PEG 8000 and the 

10 standards were calculated (using Eq 2) and four standards I4A, 14B, 14C and 14D 
were found with Proximity Index values less than the Minimal Index From the 
properties of these standards, the properties of the unknown PEG 8000 were 
calculated by averaging with (as shown \\\ Table 14 1) excellent results for all which 
were obtained in less than I minute, enabling maintenance or immediate correction of 

1 5 the level of ethylene oxide used in order to maintain the tinal quality of the products. 
Furthermore the ditVerences obtained between the results obtained by the above 
calculations and by standard methods were all inside the limits of reproducibility of 
those reference methods, namely 0.7% for the viscosity (by ASTM D44.S) and 5.8% 
and 3.6% respectively for the hydroxyl index below and above 100 (ASTM D4274). 

20 The Molecular Weight was obtained directly from the hydroxyl index with the same 
reproducibility values 

Other properties of the polyethylene oxylated butanol may be obtained in a 
similar way. 

25 : 



30 
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TABLE 10 

On line Determination ol the Properties of a High MW Polyisobutene 



In this Table 1 56E-07 means I 56x10"'^ 
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l AULt: 1 1. 1 

On line Determination of (lie Properties ol a Polyisobiiteneof Low Mol. Wt. 













1 ie 




LIS 


LIS 























ProK 

Irvlex 






4 ;(4e 07 


2 28E 07 


— n c 

.65fc-07 


9. 18E07 


1.07E-O6 


1.68E-06 


Wivelan 


nth 




















X. (f>ml 


















6079 


164S 






0 0350224 


0 034809 


0 035165 


0 Q34Q77 


0.03S02 


0,035191 


6109 


1637 




0.033505 


0 O33070C 


0,033428 


0 033629 


Q 033536 


0.033725 


0.033835 


6165 


ifff^ 


— - 


0,02964 2 


0 0296054 


0,029568 




n n j'ad7g 
0.029475 


0.029625 


0.029SS3 






— - 


0.026851 


0 02668 


0,026667 


0 0?e9S^ 


0.026525 


0.026706 


0.026544 


6215 


"~~~T609 


— - 


0.025316 


0 025199 


0.025198 




0.025055 


0.025241 


0,025087 


64 ?8 






0.020251 


0 0201638 


0.020177 


0.020325 




0,020202 


0.020O69 




1658 


— T 
— - 


0.0O894 2 


0 00905 


0 009014 


0.009098 




0.009074 


0.009016 


6532 


T53T 


— - 


0.OO791G 


0 00802G2 


0 007969 




0.008046 


0.008052 


0.008035 






— - 


0. 014038 


0 0141 226 


0 014057 


0.014102 




0.014208 


0.014125 


6696 


1493 


— - 


0 01827 


0 oi8:<04a 


0 018253 


0 01832 


0.018273 


0.01 8369 


0.01 8309 


6821 


1460 


— - 


0 023G37 


0 02:t0528 


0 023648 


0.023639 




0.O23643 


0.023666 


6901 


1449 


— - 


0 0306S1 


0 0307184 


0.03O7G7 


0.030664 


0.030755 


0.030708 


0.030698 


6925 


144T 





0 033186 


0 033I57C 


0 033188 


0.0331 15 


M 0 033189 


0.033208 


0.033088 


6964* 


1436* 


— - 


0.035052 


0 0350774 


0.0351 


0.034983 


0.035168 


0.035051 


0.03SO85 


6996 


rrtT" 


— - 


0 0345 1G 


0 0345532 


0.034541 


0 034482 


0.034606 


0.0345 


0 034637 


70S3 


j— r 




0.039698 


0 0396944 


0 039673 




0.039765 


0.039724 


0.039788 









0 041541 


0 04 1 5 304 


0 041504 


0.04134 7 


0.041604 


0.041588 


0.041609 


7092 


mTo 


— - 


0 045553 


0 04 55GUG 


0 045523 


0.045335 


0.045654 


0.045617 


0.045714 


7148 


1399 


— - 


0.052377 


0 052487 


0 052493 


0.052322 


0.052627 


0.05 24 25 


0.052568 


7158 


TsiT 


— - 


0.053934 


0 0539828 


0 054OO5 


0.053795 


0.05414B 


0.053874 


0 054092 


7199 


1389 





0.056665 


0 0569G14 


0 05701 


0.056767 


0.057037 


0.056838 


0.057155 


7220 


-rrr- 




0.056997 


0 05G99G 


0 057043 


0.056965 


0.057052 


0.056858 


0.057062 




L2|^ 




0.056741 


0 0060584 


0.05 G674 


0.056626 


0.066696 


0.056527 


0.05677 


7262 


\jT~ 


— - 


0.063595 


0 00,15540 


0 003G1 1 


0.053499 


0.053576 


0.053515 


0.053572 


7294 


L 


— - 


0 04739? 


>i ()4 7:);i9G 


0 047407 


0 047357 


0.047322 


0 047342 


0 04727 


7331 


^Tfi~ 


— - 


0.039219 


0 03914 1 


0.039240 


0.039188 


0.039132 


0 039186 


0.038954 


7348 






0 0353G9 


0 0:t52522 


0 035363 






0.035285 


0.03506 








0.02G707 


0 020097 


0.02C71 


0.026679 


0 02654 


0.02663 


0.0264 26 


2*22. 


1351 




0.O1732B 


U 01 7232 


0 017357 


0 017326 


0.017195 


0.017258 


0.01 7022 




1 266 


— - 


0.002464 


0 0024384 


0.00244 


0.002455 


0.002472 


0 002381 


0 002444 


BOOO 






0 007208 


0 007183 


0 00719 










8097 


1235 




0 017338 




0 01727G 


0.01724 


0 017281 


0017216 


0.017346 


8197 


1220 




0 035357 




0 035245 


0.035151 


0.035316 


0 034991 


0 035336 


8217 


1217 




0 0394 33 


0 0393774 


0.03944S 


0.039195 


0 03959 


0 039126 


0 039527 


8251 


1212 




0 047005 


0 04G8BHG 


0 046848 


0 046753 


0 047075 


0 046638 


0.047129 


8278 


1208 




0.05312 


0 0530G7 2 


0 053062 


0 052678 


0.053291 


O052931 


0.053174 


8333 


1200 




0 081058 


0 0R09 7 20 


0081013 


0 080977 


0 080954 


0 081039 


0.08068 


8361 






0 093478 


0 093031 


0 093501 


0 093565 


0.093469 


0 093651 


0.093469 


8375 


1194 




0.094014 


0 0941538 


0.094091 


0.094159 




0.094226 


0.094153 


8382 


1193 




0.093212 


0 09340 


0 093337 


0 093337 


0.093457 


0.093559 


0.09356 


8403 


1190 




0 0897C 


0 08992SH 


0 089B0C 


0 089808 


0.089953 


0.089966 


0 090096 


8418 


1 188 




0.08078 




0 080901 


0 086907 


0.086959 


0 087037 


0.0871 25 


8503 


1 176 




0 064 1 9 


0 004 1890 


0 064 2 1 7 


0 064345 


0 064015 


0.064341 


0.06403 
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TABLE 11. 1 (com) 
On line Detei minntiou of the Properties ol a Polyisobutene of Low Mol. Wt. 







Wt 


Meatuted 


,A 

Eviim.Metl 


IIB 


1 1C 


1 10 


1 IE 


1 IF 




Is 






<J 34E-07 


2 28E 07 


e.65e-07 


9 18e07 


1.07E-06 


leeeoe 






















A. (cm- 1 ) 




















8540 


117) 


1 


0.054912 


0 0548696 


0 054942 


0 055027 


0.OS4736 


0 05501 7 


0,064626 


esse 


1163 





0.040824 


0 040804 G 


0 04087 


0 040884 


0.040724 


0.040918 


0.040627 




1 165 














0.033613 


0 0 33462 


8703 


1149 




0.029161 


0 0291404 


0.029166 


0.029278 


0.029022 


0,029172 


O.O29074 


8726 


1 146 




0.025375 


0.025314 


0.025338 


0,02541 


0.025223 


0.025394 


0.02520S 




1136 




0.011761 


0.01 16668 


0 011764 


0.011772 


0.01 1618 


0,011658 


0,011622 




Vi.co.ny 






224,6 


221 


237 


213 


209 


243 




MN 




946 


923 


926 


930 


914 


902 


943 




LGPC 




1.59 


I 594 


1.68 


1.63 


1,55 


1.62 


1.S9 








002 


0 0216 


0.017 


0 025 


0.018 


0.021 


0.027 




TRM 




0012 


0.0112 


0.01 


0 009 




0.017 


0.01 




TRI2 




0 41 


0 407 




0 406 


0,4 25 


0,414 


0.38 




TRI2C 




0.1 12 


0 1166 


0 115 




0.121 


0,121 


0.11 




TRI2t 




0.298 


0 2914 


0.296 


0.294 


0 304 


0,293 


0.27 




TRItot 




0.578 


0576G 


0 58C 


0 557 


0.694 


0,698 


0,648 




TETRA 


1 


0 163 


0 1086 


0,169 


0 177 


0 17 


0,17 


0 167 




BUTl 
Flash 




7 0 
171 


7.3 
170 


7.2 
170 


7 2 
172 


7 5 
168 


7,13 
165 


7,47 
176 




PIBSA 




99 


99.2 


99,5 


98.6 


98 4 


99.3 


100.1 



Determination t 



TABLE 12.1 

the Properties of:! Polyethylene to Control its Manufncture 



5 3302 
4,7164 
4,2167 
3.7476 
3 2666 
2.7271 
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TABLE 12.1 (com) 

Determiiiiition of tlie Pro|)ei ties of n Folyetijyleiie to Coiitroi its Manufacture 







Wt 




2A 


12B 


12c 














Me«»uien ~ ' 

















frox. index 






0 0042G2C 


0 o67"289' 
5 


0.0095186 


0,01 281 2 


0.01 7956 


0.020326 


Wivalenoth 




















— 697 




















L£5_ 




1 386 


" r379C4 


1 3902 




1 '3076 


1 362 


1 4031 


122£. 


1 




1 2G31 


' 1 24688 


1 256 


1' 2422 


~» — i 2334 


1 2363 


2665 




1690 




1 1 3C5 


1 1301 


1 135 


1 1 2S5 


1 nei 





1 








2 017 


2 00228 


2 007 


1 996B 


1 9968 


2 0068 


2 002 


B0S2_ 


J_242_ 
























2 35 54 





2 3502 


2 3435 


2 3506 


2 3583 


2 3449 






















_^ 80BB 




















ilS!2. 


1_235_ 








3 006 


a'oi" 


3 OQo'l 


3 0777 


^ Ifllo 


— 









,284 0 2 


3 281 


3 2858 


^"2755" 


■. J 


, 


§124^ 


LiiL 




3 5B83 


3 5913a 


3 5901 


3 58B9 


3 5B27 


3 CQCH 


3 fno 


'. 2i2£. 


1121. 




:t 9312 


— TT- 


3 9335 


3 9321 


3 9259 


3 9335 


T one 


8148 


1 227 




' <X 20'i7 


4 30BBG 


4 313B 


431 1 7 


4 3029 


4 3088 


tnT~ 


8160 


1 225 




4 6904 


4 70348 


4 7083 


4 70B5 


a 6977 


4 7Q71 


4 701 3 


8172 


1224 




5 or (37 


r 074B 


ToioG 


S oar 


5 0706 


S 0702 


c 0676 


81B4 


1 222 




5 3SB4 


r ^- 770 


r 3aa5 


r 3097 


5 3754 


S 3703 


c 


8196 


1 220 




r C4r4 


r cf-CQ^ 


t\irr 




c C663 


e cfiQn 


e cfeg 








S SQ^Q 


^ f- 7 ^ 


c -) 




c; fi, 7> 


1 B 


5.5976 


8304 


lio4 




62, 


• aii'^'^ 


77r 


, ^ 

3 9785 


Qfi^ 

3 9664 


qI,^ 

3.9677 


^^^*^ 




1203 




^ 7 97 

3 097 


3 71 178 


3 7156 


3 71 26 


3 7071 


^'Oei 


3.7145 


8328 


1201 






3 44774 


3 4551 


iAi'ie 


3 4421 


3 445 


3.4509 


8340 






3^)655 


3.185)2 


3 1931 










8362 


t 197 




2 94 26 


2 9428G 


2.9445 


^ 


7 Q ?^ 


qo^'^ 


^ nf^g 


8364 


11 96 




2.7439 


2 74022 


2.7361 


2 7369 


2 7463 


^2 727 


2 7548 


8376 


1 194 




2 C033 


2 G0352 


2 5999 


2.603 


2 G104 


2.5913 


2.613 


8388 


1192 




2.54G5 


2.55012 


2 5457 


2 653 


2.5529 


2,5426 


2.5564 


8400 


1 190 




2.54 


2 53978 


2 5322 


2.6348 


2 6421 


2 5409 


2 5489 


8412 






2.4920 




2 4825 


2.4746 


2 4951 


2 497 


2.4991 


84 24 


1187 




2 3G:i5 


J 3577C 


2 3517 


2 :i43 


2 3666 


2 3639 


2.3636 


8436 


1 185 








2 1797 


2 1783 


2 2041 


2 1911 


2.1953 


8446 






2 047 




2 0245 


2 0367 


2 0614 


2 0337 


2 044 


84 60 


1 I 82 




1 9027 






1 B99C 


1 9027 


1 8872 


1 .8929 


8472 






1.7CC2 


1 7G79G 


1 75G2 


1 7616 


1 7656 


1.7567 


1 7509 


8484 


1 1 79 






1 00820 1 1 GGB5 


1 0062 


1 6732 


1 6694 


1 662 


8496 


1177 




1 C317 


1 G1802 1 I 0109 


1 6228 


I G257 


1 618 


1 6157 










4.3 4 2 








4.6 








0.953 


0 902 1 0 902 


0 953 


0 902 


0.952 


0 951 
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lAULL 13.1 

Detcrmiiiutioii of I'crceiiiaijcs of Products Icjwiiig the Reactor 







Wl 


13A 


138 


13C 


130 




















Prox Index 






Emmated ^^^^^ 

: ] 


0.074338 


on 208 


0.16551 


Wavel«ngth 


































0.0046638 


0 9506 




0,9167 


0.8792 


0 9062 


6620 




0,01381 1 


1 6348 


vT 


1,5997 


1.5768 


1 5994 


6652 


1503 


0 010524 


1 7996 


77 


1.784^ 




1 7742 






0,015252 


2 2127 




2.1896 


2 1635 


2 1819 












2.2957 










0,0051909 






2.7882 










0,0024662 






















3.7794 






















2226 






3 9273 














0 14G04 




— 


4 1628 


"4" 1 702 


4 2177 


















1263 








_ 










136; 


■ -Q 33437 


' 6563 




1 7027 


1 5795 


1 598 


7466 


1340 


0,025072 


° ^ 




0 4291 


0 3898 


0 3968 








2_Zi. 


^ 

— 


0 7664 


0 654 2 


0 6631 


8100 




0 79I37 






3 0876 


3 1733 


3 2087 


8260 


121 2 


1 5868 


0 51 07 




G 4577 


6 689 


6 7332 


8332 


1200 


0 99957 






6.9601 


7 1546 


7 1455 


8434 


iiec 


0 62305 


G 3714 




S.S339 


0.3954 


6 3456 


8592 


1 1C4 


1.9117 


4 31:11 




4 4163 


4 2975 


4 2704 


8660 


1 165 


0 G550i7 


:< 0575 




3.7572 


3.7107 


3 7241 


8710 




0 594Cr. 


G 0O15 




5 0538 


5 0873 


5 0805 


8767 




0 01 289 






5 3749 


6 5012 


5 4754 


8796 




0 70C3B 


0 598B 




5.5164 


5.6277 


5.5518 


8815 


1 1 34 


0 56808 


4. 0700 




4.5946 


4 6884 


4.6 


8941 


1 1 31 


0.34646 


J 4155 




3.3726 


3.4198 


3.3602 


8860 


1 1 29 


0,25730 


2.8157 




2.7751 


2 8008 


2.7556 


8936 


1119 


0 19092 






1.2403 


1 178B 


1 1695 


8955 


1117 


0.174 34 


0 9454 




0,9796 


0.9423 


0.9195 
























12.85 




16 99 


12.24 


12.39 


Heavy 






42.71 




3G.19 


41.94 


40.78 
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TABLE 13.2 
Nov "Staiulard" MCI obl:iiiu*ii by Deusiricntion 



0 B792 

1 576B 
1 7564 
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TABLE 13.3 
New "Standard" MC2 obtained by Deusificatioit 







MC2 


13B 


130 


13E 


Fraction in the 

Miilurc 






0.877 


0.136 


0.259 


WaveleiMith 










k (cm-1 1 


iiii.ii 










C4 60 




0 93676 


0.9167 


0 90C2 


0.9887 


6620 


161 1 


1 6276 


1.5997 


1 5894 


1.7026 


G6B2 


1503 


1.8053 


1 7841 


1 774 2 


1.8607 


671 1 


1490 


2.2059 


2 1896 


2 1819 


2.2483 


6730 


1486 


2 3126 


2.2957 


2 2865 


2 356 


6790 


1471 


2 796 


2 7882 


2 7702 


2 8087 


GB24 


1460 


2.8976 


2 8978 


2 seoG 


2,8877 


6996 


I4 2H 


:» 784 2 


3 7794 


3 8055 


3 8116 


7028 


14 2:( 


4 Bbb2 


4 ei4G 


4 9319 


5 1604 


707G 


14t:t 


3 9406 


3 9 644 


3,9604 


3 8653 


7150 


1399 


4,1696 


4 1628 


4 21 77 


4,2177 


7215 


13H0 


2.6017 


2 6289 


2,638 


2.5287 


7263 


1377 


3 1154 


3 1596 


306781 2 9409 


7344 


1362 


1 6437 


» 7027 


1 598 


1 4198 


7465 


1340 


0 4 3413 


0 4291 


0 3988 


0 4326 


7504 




0 70472 


0 7664 


0 6631 


0.474 


aioo 


1235 


J 0874 


3,0875 


3 2087 


3 1509 


8250 


1212 


e 5187 


6 4677 


G 7332 


6,8378 


8332 


1200 


7.0364 


6.9601 


7 1455 


7,3483 


6434 


1 IBO 


642 


6.5339 


6 345G 


5 9953 


8592 


iir.4 


4 3145 


4 4163 


4 2704 


3 9466 


86G0 




3.0556 


3,7672 


3 7241 


3,3484 


8710 




4.9937 


6 0538 


5 0805 


4 8367 


87G7 






5 3749 




5,5484 


8796 


1137 


5 5797 


5,5154 


5 5518 


5,7792 


8815 


1 134 


4 C687 


4 6946 


4,6 


4,8836 


8841 


1131 


3 4 355 


3 372G 


3.3602 


3.6088 


8860 


1129 


2 84 


2 7751 


2.755G 


3.0153 


89 3C 




1 2142 


1 ,2403 


1,1695 


1 1023 


8956 




0 9661 


0 9796 


0 9195 


0 8959 














L.nt.t Fr.,ci,on 




13 20 


16 99 


1 2 39 


0 17 






4 1 99 


3C 19 ! 40 78 


CO 98 
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TABLt 13.4 

Determination of Percentages of Product leaving Reactor 
based on New "Standards" 



W«vl«0(Hh 
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TABLE 13.5 

Iiiimeiiiiite use o( iu-m "Siniuhii J" by AutOiidjiistinent 



W»vel«ntith 



0 004GG3B 
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lAULE 14.1 

Determination of Properties ui a PEG to control u Production Unit 







_W|_ 




li 


— mi — 


111! 




14E 








Mcawircl 


r^limaUAl 










uu ^ 


r»dex 






0.O0OO6519 


0.00000977 


0.00007494 


0,00010457 


0.0001 1674 


jvej 


tnoth 


















*• 
















4164 


2122. 




Q^^^^^ 


o'^^yGO^ 










0.11846 
















0.23696 


0.23729 


1221. 


2210 




0.054854 


0.O5G7O4G 


0 054802 


0.056171 


0.057599 


0.058246 


4636 


2068 




0 048031 


0 0465855 


0 048451 


0.048238 


0.048588 


0.049065 




1933 




0 059102 


0.06188575 


0.061585 


0.061986 


0.064102 


0.069B7 


5436 


1840 




O.0S6O88 


0.0555625 


0.056049 


0.0555 


0.055054 


0.056647 


5544 


1804 




0.068293 


0.06745175 


0OC7773 


0066995 


0.067073 


0.067966 


574B 


















58S6 


1708 




0.037196 


0.0374 250 


0 037003 


0.037G03 


0.037295 


0.037BO1 


6624 


1510 




0014687 


0 015498 


0 014728 


0 01C063 


0.01551 


0.015691 


6664 


1496 




0 019234 


0 01960475 


0 01 9008 


0019792 


0.0I9B26 


0.019733 


6720 


1488 




0.020139 


0 02084475 


0 020349 


0.021013 


0.02094 2 


0 021075 


6792 


1472 




0.022761 


0 02art5:t 




0 023C3C 


0.023768 


0 02357 


6972 


1434 




0 030389 


0 0307H3.'5 


0 0301 75 


0.030929 


0.031 19 


0 030839 


7092 


1410 




0 021 148 


0 02100025 


0021 221 


0.022095 


0 021621 


0.021 724 


711G 






0 021822 


0.0224GI 75 


0 02158 


0.022824 


0.O2273 


0.O22713 


7920 


1263 




0 0O4193 


0 005220:i5 


0 004 21 GG 


0 0057718 


0 0OS4551 


0.0054389 


8172 


1224 




0 03318 


0.0330O15 


0 0a2O74 


0 032927 


0.032945 


O.034O6 


8352 


1197 




0 011492 


0.01 1tt575 


0 01092G 


0012427 


0.01 1449 


0.012628 




















HydroKvl Index 




13.1 




13 1 


13 


13.2 


13.16 






701 


70.1 


703 


710 


708 


696 


Molec. Waiaht 




S560 


KG12.5 


H400 


8370 


8450 


9230 
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Examples 15 

Determination of the Properties ot a base Oil 

The NIR spectaim of a base oil D which was a 500 neutral oil was measured 

between 4800 and 4000cm-l with normalisation of the absorbances, [the base line 
5 being taken at 4780cm- ' ] 

NIR spectra were recorded for a series of standard base oils of krwwn properties. By 

the Minimal Index Procedure described above, with use of Equation 2, and non 

weighting of the absorbancies. the Minimal Index was calculated to be 5x10"^. 

Following reference to the bank of data on the standard base oils, 3 standards 15A, 
10 15B and 15C were found with the proximity index with respect to oil D less than 

5x10-7 The propenies of those standards and their spectra are shown in Table 15.1. 

By averaging the properties of the standard samples, various propenies were obtained 

for the oil D The Table 1 5 shows the estimated properties as well as measured 

properties of oil D for comparison. 
15 All the properties were determined in a single analysis and without any 

regressional type calculations and wiih a degree of precision in line with the 

reproducibility of the reference metlHids 

Other properties ofD can be determined in a similar way 

Example 16 

20 Determination of the properties of a Process Oil 

The method of Example 15 was repeated witli a process Oil Reference I6D of 
"Enerthene" type which was a mixture of neutral base oil and aromatic-containing 
vacuum distillate extract From a bank of standard process oils of this type. Minimal 
Index was found by the Minimal Index Procedure to be 5x1 0'^. providing a sphere of 

25 identity. 3 standard oils 16A. 16B and 16C were found with proximity indices with 
respect to 16D less than the Minimal Index and hence inside that sphere The 
properties of oils 16 A, 16B, I6C and their spectra and the spectrum of 16D are given 
in Table 16.1. By arithmetic mean averaging of the properties of 16 A, 16B and 16C, 
the properties of I6D were estimated, and these together with the measured properties 

30 of 1 6D are given also in Table 1 6 1 

The single analysis gave all the propeities without regression calculation and 
with an accuracy in line with the leprtHlucihiiiiy of tiie leference methods. Other 
properties can be determined in a .similar way 
' In Table 16 I, the expression 4.20 E-04 means 4 2x10"' and PC A means 

35 Polycyclic Aromatic hydrocarbon 
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Example 17 

Determination of the Pronenies of i\ C rude Haratfin 

The method of Example 15 was repeated with a crude paraffin (a "slack" wax), 
which was a mixture of paratiin wax and base oil called Gatsche D From a bank of 
5 standard crude paratfms and their propenies and spectra, the Minimal Index was found 
to be 5x10-5, by the Minimal Index Procedure. 3 standard crude paraffins Gatsche 
17 A, 17B and 17C were found with proximity indices with respect to I7D inside the 
sphere of identity The propenies of paratiins 17A. 17B. 1 7C and their spectra, and 
the spectrum of 1 7D are given in Table 17 1 By arithmetic means averaging of the 

1 0 propenies of 1 7A, 1 7B and 1 7C the propenies of 1 7D were estimated, and these 
together with the measured propenies of 1 7D are given also in Table 17 1 The 
properties determined were density, viscosity at lOOX and oil content. The 
determination of the oil content of crude paratfin according to the above procedure is 
particularly remarkable because no on line process capable of measuring this is known 

1 5 today 

All the properties were determined in a single analysis without any regression 
type calculation and with an accuracy in line with the reproducibility of the reference 
methods. Other propenies of D can be determined in a similar way 
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TABLE 15.1 
Detenniimtion oT the nroDerties of Base Oil 



Olio 



OilD 
Esliitial 



X(an-1) X(nm) 



5.0333333E-O4 



S.OOOOOOOE-04 



5.70O0O0OE-O4 



5.80000006-04 



5.9OO00OOE-04 



5.9000000E-04 



6.1000000E-04 



6.3000000E-04 



6.1O0O000E-O4 



7 80000006-04 



7.700OO0OE-04 



7.60000006 -04 



84000000E04 



6 40000006 -04 



8.3000000E-04 



8 466€6€7E-04 



8 4000000E-04 



e.5000000E-04 



3.10OO00OE-O4 



8 0333333E-04 



8.0000000E-04 



B.1000000E-04 



8.6000000E-04 



8 73333336-04 



8 8000000E-04 



9 3000000E-04 



9.233333364)4 i 9 2000000E-04 



9.2000000E-04 



2.00000006-03 



2.0000000E-03 



2,0OOO0OOE-03 



4.23000006-03 



5.52333336-03 



5.5300000E-03 



S.5500000E-03 



7 43O00O0E-03 



7.43OOOO0E-03 



7.4SOOOO0E-03 



4330 2309 



4240 235fi 



9,8000000E-03 



9.820O000E-03 



2.04S0000E-02 
2.6140OO0E-O2 



3.2330OO0E-02 



3.6C0O0O0E-02 
2.997OQO0E-02 
2.6S)OO00E-02 
2 454000QE 02 



2.04800006-02 



3.2136667E 02 j 3.23000006-0 



2.6130000E-02 
3.2330000E-02 



2.5390OOOE-O2 
2.914O00OE-02 
3 73800006-02 
3 2840000E-02 



3.6566667E-02 
2.9933333E-02 
2 6S66G67E-02 
2 4500000E-02 
2,383000QE-02 
2S346667E-02 
2.9163333E-02 



3.6550000E-02 
2.9930OOOE-02 
2.6570OOOE-02 
2.4490000E-02 
2.3820000E-02 



2.919OOOOE-02 
3.7330OOOE-02 



2.5a20OO0E-O2 



- 2.5860000E-02 



2.6083333E-02 
2.585O0OOE-02 
2. 56533336-02 
2.5853333E-02 



2.80a0OOOE-02 
2.6120OO0E.02 



2.6S9OOOOE-02 
2.4S10000E-02 
2.3870000E-02 
2.53400006-02 
2.91400006-02 
3.73700006-02 
3.28200006-02 
2.61000006 -02 
2.6050O00E-02 



2. 58300006-02 



2.5820000E-02 



2.5620000E-02 
2.58400006-02 
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TAliLE 15.1 (coiic) 
Deteriniiuitioii of the iirooerties of Base Oil 



Proximity Index 



Cafbon % 
lnllammat>ilil< 
Pcnsky-Manens 

JC 

Nttrogan base 



2.5820000E-0Z 
2.S970000E02 
;.62;000OE-O2 
;.66O0000£-Q2 
2.6370O0OE-02 
2 .57100006 -02 
2.46200006-02 
2.399000CE-OZ 



OilD 
EsiimaieO 



2.S646667E-02 



2.62366676-02 



2,6390OO0E-02 



2.46333336-02 



2.26OOO0OE-O2 
2.2730O00E-O2 



2.36400006 02 



2.26133336-02 
2.27400006-02 
2.33633336-02 
2 36366676-02 



2.58200006-02 
2.60100006-02 



2,62400006-02 



2.63800006-02 



2.59700006-02 
2.61900006-02 
2.6S80000E-02 
2.63700006-02 



2.3980000E-02 
2.30500006-02 
2.26000006-02 
2.2720000E-02 
2.33300006-02 
2.3650000E-02 



2,66400006-02 



2.5730000E-O2 



2.4640000E-02 
2.40000006-02 
2.30600006-02 
2.26100006-02 
2.2730000E-02 



2.36100006-02 



2.40000006-02 
2.30900006-02 
2.26300006-02 
2.2770000E-02 
2.33900006-02 
2.36500006-02 
2.07400006-02 
1.5190OO0E-02 
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I ABLE 16. 1 
1 ilic t)ni|>ci tics »r ;i im occnn oil 



Olio 



3.S0OOOOOE-(M 



3.500O000E-O4 



7.300OO0OE-O4 



6.466666 7E -04 



6.6000000E-04 



9.60CXJOOOE-04 



9.4000000E-04 



9.9000000E.O4 



2.160OOOOE-O3 



2.09OO0OOE-03 



2.0900000E-03 



2.1?OOOOOe-03 



2.51000OOE-O3 



2.43000OOE-03 



2.41000006-03 



2.47000006-03 



2.9500000E-03 



2.8600000E-03 



2.9100000E-O3 



3.00000006 -03 



3.110000OE-03 



3.0966667E-03 



3.0100000E-03 



3.080O0OOE-O3 



3.20O0O0OE-O3 



3.1833333E-03 



3.1200000E-03 



1.1 8000006.03 



3.2500000E-O3 



4258 234S 



3 0700000E-03 
2,6300000E-03 



2.3200000E-03 



2 3400000E-03 



6.140O00OE-O3 



2.0970000E-02 



3.3140000E-0 



2.85800006-0 



2.S410000E-02 
2. 39300006-02 



3.1466667E03 
3.0S0OOOOE-O3 
2.SBOOOOOE-03 
2.3133333E-03 
2 1933333E03 



2.3133333E-03 
3 20OOOOOE-03 



7.6100000E-03 



1.9386667E-02 
2.0963333E-02 
2.19n333E-02 
2.2530OO0E.02 
2.3053333E-02 
28183333E-02 
3.3196667E-02 
3.87aOOOOE-02 



2.9810000E-0: 



2.5393333E02 



3.1100000E-03 
3.OO0OOO0E-03 
2.SS00O0OE-O3 
2.2S00000E-03 
2 15Q0000E-03 
2.0O0O00OE-03 
2.29OO000E-03 
3.1SOO0O0E-O3 
6.0800000E-03 
7.780O0OOE-O3 



3.160OOOOE-03 
3.05OO0O0E-O3 



2.580OOO0E-03 



2.160O0O0E-O3 



2.29OO0OOE-03 



6.0600000E-03 



3.1000000E-03 
2.6100000E-03 
2.3eOOOOOE-03 
2.2700000E-03 
2.01OOOO0E-O3 
2.36000006-03 
3.2600000E-03 
6.1500000E-03 
7.88000006-03 



2.187O00OE02 
2.2510000E-02 



2.3030000E-02 
2.8170000E-02 



2.0940000E-02 
2.1930000E-02 



2.2520000E-02 
2.3030000E-02 



2.1020000E-02 
2.19400006-02 



2.2560000E-02 
2.3100000E-02 



2.82200006-02 



3.8aS0000E-02 
34360000E02 
3O83O0OOE-02 
2 85600006 -02 
2,63100006-02 
2.S200OOO6-O2 
2.S8O0O00E-O2 
2.8300000E-02 
3.250OOOOE-O2 
2 9370000E-02 
2.7120000E02 



2 4000000E-02 



3.3220000E-02 
3.88100OOE-02 
3.4300000E-02 
3.08700OOE-02 
2B5600OO6-O2 
2.640OOOOE-O2 



3.e680000E-02 
3.43000006-02 



2 57eOOOOE-02 
2.8290000E.02 



2.86100006-02 
2.64200006-02 
2.52400006-02 
2.5780000E-02 
2 8200000E 02 
3.2260000E-02 



2.7170000E-02 



2.402OOOOE-O2 



2.5360OO0E-02 
2.39B0OOOE-02 
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TABLE 16.1 (tout) 
Detenniiiatioii of the properties of a process oil 







mi 


CHID 


_OiM6^_ 







Proximil tiidex 






roioo'^ — 


2±2S. 


-LI*ll° 


°" ^ 

i:12lI2 
















).(cm-1)^(nm) 














"^^ 




2.2630000E-02 


2,26300006-02 


2,26800006-02 


2,26500006-02 


2.25600006-02 


4200 2380 




2.20€0OO0E*O2 


2.2106667E.02 


2.2140000E-O2 


2. 2 1600006 -02 


2-20200006-02 


*^^^ 




2.2010000E-O2 


2.2043333E-02 


2.21100006-02 


2,20700006-02 


2.19500006-02 


— -"- 




2.2220OO0E -02 


2.22266676-02 


2,23100006 -02 


2,22300006-02 


2.21400006-02 


_4176 — 2394 





2.2780000E-02 


2.2816667E-02 


2,28600006 '02 


2.2840000E-02 


2.27S0000E-02 


ills — 





2.3160000E -02 


2.321 3333E-02 


2.3290000E-02 


2,32100006-02 


2.31400006-02 


4160 2403 





2. 284 00006-02 


2 2SS0OOOE-02 


2.28900006-02 


2,28600006-02 


2.2800000E-02 


^""^ 


-1 


2 18100006-02 


2 1S43333E-02 


2.19O00OOE-02 


2. 18600006-02 


2,17700006-02 






2.0630000E-02 


2.O630OOOE-O2 


2.0700000E-02 


2,06400006-02 


2,05500006-02 


^^^^ 




2.O1700OOE-02 


2.01866676-02 


2.02400006-02 


2.02200006-02 


2.01000006-02 


4104 2436 


1 


1 9520O00E-02 


1.95633336-02 


1.95900006-02 


1,95900006-02 


1,95100006-02 


4092 2443 




1. 95300006 -02 


1. 95933336-02 


1.9640000E-02 


1.9610000E-02 


1 .953O00OE-02 


4080 2450 




2.1540000E-02 


2 15133336-02 


2,15500006-02 


2.15300006-02 


2,14600006-02 






2 .3530000E-02 


2.35300006-02 


2.35500006-02 


2.35300006-02 


2.35100006-02 


406« 2458 




2 34300006-02 


2 34433336-02 


2,34600006-02 


2,34500006-02 


2,34200006-02 


4048 2470 




1.89900006.02 


1.9O10O0OE-O2 


2.90500006-02 


1,90200006-02 


1.89600006-02 


4000 2500 




1.46300006-02 


1.4593333E-02 


1 ,46200006-02 


1.45800006-02 


1.4S80000E-02 
















0«nsity kg/I 




0.9348 


0.9351 


0.9350 


0.9345 


0.9360 


Sulphur % 




2.25 


2.28 


2.37 


1.98 


2.51 


PCA% 




2.60 


2.57 


2.88 


2,74 


2.1 


100-C sCl 




33.19 


33.22 


33.23 


32,25 


34.18 


Flash Poinl 
Clevaland *C 




310 


311 


310 


308 


315 



wo 96/1 1399 




PCT/GB9S/023S3 



lAULt: 17.1 

Detenninatioii ol (lie properties of a crude parafTiii 



Ga(tche 0 
Etiitiiaied 



Gatsche 17A 



7.0000000EOS 



a.OOOOOOOE-M 



a.oooooooEos 



8,00O0O00E-O5 



1,9000000E-0< 



2.0O00OOOE-O4 



3.0000000E-(M 



7.6666667E-0S 



7.0000000E-0S 



2.7OOCX)O0E-O4 



V9000000E-04 



I.90OOOO0E04 



3.0333333E-0 



2.4000O0Oe-O4 



2,600O0O0E-O4 



2.0333333E-04 



3.4333333E-04 



3 6000000E-04 



3 000OO00E-O4 



3.7000000E-O4 



5.30000006-04 



3.9333333E-04 



3.80O0O0OE-04 



4.0000000E-04 



7.70000O0E-04 



3.0000000E-04 



3 OOOOOOOE-04 



2.9O00O0OE-O4 



5.5000000E-04 



2.S000O0OE-O4 



2.6000000E-04 



6.400OOD0E-04 



2. 90000006 -04 



3.0000000E-04 



6.1000000E-04 



S,000O0O0E-O4 



B.9000000E-04 



7.3666667E-04 



7.2000000E-04 



4440 2252 



3 0766667E-03 



3.1900000E-03 



30400000E-03 



3.0000000E-03 



.1400OO0E-03 



s.saooocoE- 



5.5S00OOOE-O3 



5.S50OOOOE-03 



5.3400000E-03 



4416 2264 



7.4700000E-03 



7.3900000E-03 



7.6000000E-03 



7.4200O00E-03 



9.8500000E-03 



4392 2275 



1,720OOO0E-O2 



I.64220OOE-02 



2.7510OOOE.02 



2.7420000E-02 



3.7070000E-02 



3.7120OOOE-02 



3.7170000E-02 



3.38S0QO0E 



3.391OOOOE-02 



3.4130OO0E-02 



34660000E02 



2.5860000E-02 



2.S900000E-02 



2.S9B0OO0E-02 



2.2440000E 



2.2710000E-02 



2.31OOOO0E-O2 



2inOOOOE-02 



2.126OOO0E-02 



2.126O0OOE-02 



2.1310OO0E-02 



4280 2336 



4248 2354 



2 5370000E-02 



2.6273333E-Q2 
2 S383333E 02 



2 4510000E 02 
3.17O0OO0E-02 



3. 2 7600006 -02 



2.S330000E-02 



2.4760000E-02 



2.6260000E-02 



3.2960000E -02 
2.65DOOOOE02 



2.S410OOOE02 
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TABLi: 17.1 (com) 
Determination ul ilie properties ol'n crude parafTtn 



Galsche D 
MeasuiiKl 



2.69«10-^ 



2.6390000E-02 



2.6<93333E-02 



2.6390000E-02 



2.68S0000E-O2 



2.6a9000OE-O2 



2.784OOO0E-02 



2.7633333E-02 



2. 76100006 -02 



2,7640000E^2 



2.7650000E-02 



2.81900006 -02 



2.8150000e-02 



2.8210O00E-O2 



2,8140O00E-O2 



2.821QOOOE-02 



2.8200000E-O2 



2.854OO0OE-02 



2.8663333E-02 



2.873OOOOE-02 



2.86000006 02 



2.8970OOOE-02 



2 9O00OO0E-O2 



2.9010000E-02 



2.9030000E-02 



2.841OOOOE-02 



2.8526667E-02 



2.8S50000E-02 



2.767OOOO€-02 



2,7740OO0E-O2 



2.7670000E02 



2.7760000E-02 



2.679OOO0E-02 



2 6a40000E-02 



2.6660O00E-O2 



2.685OC00E-O2 



2.57200006-02 



2.56600006-02 



2.5660000E-02 



2.S660000E-02 



2.457OOOOE-02 



2 4S06667E-02 



245100006-02 



2.4420000E-02 



2.4590000E-02 



2.40100006-02 



2 39533336-02 



2 39700006-02 



2 39100006-02 



2.39800006-02 



2.392O0OOE-O2 



2.40566676-02 



2.39700006-02 



2.40700006-02 



2.4480000E-02 



2.45aOO0OE-O2 



2.441000OE-02 



4068 2458 



2.45300006-02 



2.4480000E-02 



2 .45600006 -02 



2.0460000E-02 



2.0390000E-02 



2.048O00OE-02 



2.0240000E-02 



2.04S0OOOE-O2 



4000 2500 



1 .S6600Q0E-02 



Oil Content K 
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Claims 

I A method of determining or predicting a value which is a value ofa 
property ofa material X or a propeny ofa product ofa process from said material or 
yield of said process, which method comprises measuring the absorption Djx of said 
material at more than one wavelength in the region 000-2600 nm, comparing said 
5 signals indicative of said absorptions or a mathematical function thereof with signals 
indicative of absorptions Djni ai the same wavelengths or a mathematical function 
thereof for a number of standards S in a bank for which the said property or yield P is 
known, and choosing from the bank ai least one standard Sp^ with property P,^^ said 
standard having the smallest averagL- value ol the absolute ditference at each 
1 0 wavelength i between the signal for the material and the signal for the standard S,-n to 
obtain Fx. with averaging of said propei iies or yields P^ when more than one standard 
Sm is chosen 

2. A method according to claim I comprising comparing absorptions Djx (or a 
derivative thereof) with absorption Djiii or a derivative thereof 
15 Ti. A method according to claim i or 2 wherein the standard Sj^^ chosen for the 
property or yield wanted is such that in relation to the unknown material X and each 
chosen standard the following function is met 
'xm experimental error 

20 IDix 

wherein ixm is the proximity index and is deliiied by \-(\m) = KDiv;-D|,-„)2 and the 
experimental error is in determining said piopeny oi yield m the standard 
4, A method according to claim wherein the proximity index is less than the 
minimal index i,,i which has been dcici ninied front the standards Sa. Sb. Sc by (a) 
25 calculating for each pair of standards Sa/Sb. Sa/Sc the value of i^Ca.b) etc. (b) relating 
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the values of i-(a,b) etc to the corresponding differences EP (ab) in properties Pa. Pb 
etc (c) calculating tor each value L for which i2(ab) is <L, the average of the 
corresponding differences EPab, (d) calculating Minimal index from the value of 
minimal index i2(ab) where average EPab is the same as reproducibility standard for 
5 the propeny 

5 A method according to claim I, wherein the properties of synthetic standards, 
which are mixtures, and their spectra for consideration for possible choice for are 
estimated from existing standards in the bank for which, in respect of each existing 
standard for use in said mixture equation (4) and (5) are met, 

1 0 MinCj - uj < Cij < Max Cj + uj (4) 

and I Cij = 1 (5) 
wherein Cjj is fraction of component j m the sample i 

Min Cj is the minimum of j in the standards in the bank or in the samples for which the 
method is to be used 

1 5 Max Cj is the maximum of j in the standards in the bank or in the samples for which the 
method is to be used and uj is between I U and 0 05 

6 A method according to claim 4 wherein properties of Standards and spectra for 
consideration for possible choice are estimated by interpolation from measured 
properties of Standards and spectra for whicli the proximity index with respect to the 

20 unknown X is not more tlian 10 times the Minimal IndeN 

7 A method according to claim 5 wiierein at least one of (i) the estimated 
Standards and the corresponding specu a. and (ii) tlie propeny l\ of the unknown 
material and its spectrum, are added to the biink 

8 A method according to any one of tiie preceding claims wherein the property 
25 Pj^ or yield is compared to the desired value and any deviations used in a closed loop 

control system to control the processing equipment in relation to a process for which 
the material is a product or a feed 

9. A method according to any one of claims 1 -8 wherein the property is a 
physicochemical propeny of material X 
30 10. A method according to any one of claims 1-8 wherein the property is a 
physicochemical propeny or yield of a product of a process to which at least one 
material X is a feed. 

11 A method for adding an extra synthetic standard to a bank of known standards, 
each of which relates at least one absorption in the 000-2600nm region (or signal 
35 indicative thereof or mathematical function of said absorption) of a known material to a 
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known propeny related to that material, which propeny is of said material which is a 
feed to a process, or is of a product from said process or yield of said process, which 
method comprises choosing from the bank at least 2 standards for which equations 4 
and 5 as defined in claim 5 are met. considering mixing the chosen standards in at least 
5 one proponion to produce at least one mixture for use in a synthetic standard, and 
estimating the spectaim of said mixture according to equation 6. 

SMi = I CijXSj (6) 
where Sj is the spectrum in the mixture of component j in the calibration matrix, 
and estimating a propeny of said mixture according to equation 7 
10 PM.= I(^ijXP| (7) 

where Pj is the property of componeni j. 

and then adding; the specmim and pi upcriv of each "n^ixlure" to the bank . and using 
them in at least one model iiivolvin\ji ;i correluiiou/regresision approach to relate NIR 
spectra to at least one property 
15 12. A method according to any one of the preceding claims wherein said process is 
a hydrocarbon conversion or separation process, preferably a reforming or catalytic 
cracking or hydroireatmenl. oi' distillation oi blending. 

13 A method according to claim 12 wherein said process is a distillation to give at 
least 1 distillation product and a residue and the properties/yields are obtained in 

20 respect of said product and/or residue 

14 A method according to any one uf the preceding claims wherein said property 
is in respect of a motor fuel and is at least one of an Octane Number, vapour pressure, 
volatility percentage distilled at 70 and at 100°C, gum content in mg/lOOmi and 
content of sulphur, benzene oi nieihyl len butyl ether 

25 15. A method according to claini 14 wherein the propeny is in respect of a blend 
comprising gasoline, the spectra are nieasuied on feeds to said blending, and by 
calculation the blend index tibiained as a hnear or non linear function 

16. A method according to any one of claims 1-13 wherein said propeny is in 
respect of gas oil and is at least one of cetane index, cetane number, percentage of 

30 sulphur, density at \^°C, clear pomt. cloud point, filterability and viscosity at 40°C. 

17. A method according to any one of claims 1-11 wherein said propeny is in 
respect of a crude oil and is at least one of density, percentage of sulphur, viscosity at 
100°C, content of paraffin and residual carbon percentage (Conradson Carbon) 

'18 A method according to claim 12 wherein said propeny is in respect of a feed to 
35 a reforming process and is at least one of percentages of saturated linear, isoparaftins. 
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napthenics, and aromaiics and density 

19. A method according to claim 12 wherein said propeny is in respect of a feed to 
a fluid catalytic cracking unit and is at least one of the density, the weight percentage 
of sulphur, the aniline point, viscosity at 100°C, refractive index at 20°C or 60''C, 50% 

5 distillation point, molecular weight, percentage of aromatic carbon and the KUOP, 
crackability or cokability of the feed or yield of gas, gasoline, gas oil or residue. 

20. A method according to claim 12 wherein said percentage is in respect of the 
feed to a hydrogenation unit and is at least one of percentages of linear saturation, 
isoparaffins, napthenes, linear olefms, cyclic olefins, benzene, toluene, xylene. 

10 alkylbenzene. density, or yield of light cut, heavy cut, or raffinate or benzene. 

21. A method according to any one of claims 1-11 wherein said process is at least 
one of a polymerization, an oligomerization or an organic reaction in which at least one 
of the reactant and a product is a functionalized compound 

22. A method according to claim 21 wherein said process is a polymerization. 

1 5 23. A method according to claim 22 which is an addition polymerization of at least 
one unsaturated hydrocarbon, and the property is at least one of number and weight 
average molecular weight, molecular weight distribution, viscosity, viscosity index, 
fluidity index, density, chemical composition such as percentage of at least one 
monomer, or unsaturation or side chain group, crystallinity, rigidity, a flow property . 

20 draw strength at the Mow threshold, cracking resistance and shock resistance 

24. A method according to claim 2."^, wherein the process is a polymerization of at 
least one alpha olefin of 2-8 carbons and the property is of the product and is at least 
one of the density, tluidity index, degree of conversion, content of volatiles and. in the 
case of copolymerizaiion percentage of comonomer. 

25 25. . A method according to claim 2.'> wherein the process is a polymerization of 

isobutene and the propeny is of the product and is at least one of the viscosity, number 
average molecular weight, distribution of molecular weights, inflammability point, 
content of butene 1 unsaturation and maleinization index, and percentage of light and 
heavy fractions from distillation of the direct product of the process 

30 26. A method according to claim 22 wherein the process is the polymerization 

reaction of an epoxide in the presence of an organic compound containing at least one 
hydroxy group and the propeny is at least one of the degree of conversion, and of the 
product, hydroxyl index, viscosity and molecular weight 

27. A method-according to claim 2 1 wherein the proce.ss is a reaction in which at 
35 least one of a rea'ctant and product is a functionalized compound, and is a hydration. 
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dehydration, etherification, esterification, oxidation, ammoxidation or carbonylation. 
28. A method according to any one of claims 1-10 wherein said material X is a 
composition comprising part of a lubricating oil fraction obtainable from a vacuum 
distillation of oil 

5 29 A method according to claim 28 wherein said material comprises at least one of 
a lube base oil. bright stock, process oil. wax and paraftin 

30. A method according to claim 28 or 29, wherein the property is of the material, 
which is (i) a base oil and the propeny is at least one of the density, viscosity. Viscosity 
Index. Flash Point. Pour Point, and the content of sulphur, nitrogen base, aromatic 
1 0 carbon and Polycyclic aromatic hydrocarbon, (ii) a wax or paratVm or mixture thereof 
with a base oil and the property is at least one of the density, viscosity and base oil 
content or (iii) a process oil and the property is at least one of the density, clear point, 
viscosity and content of sulphur or polycyclic aromatic hydrocarbon 

31 A method according to claim 30 wherein the propeiiy is of a finished 
1 5 lubricating oil which comprises at least one non hydrocarbon additive. 

32 A method according to any one of the preceding claims which is computer 
implemented 

33. A computer implemented method according to claim 32 involving a 
spectrometer linked to a process line containing a material X, a computer linked to the 

20 spectrometer, and a controller linked to the computer and the process line, the 

computer including databanks having stored tiierein signals indicative of absorptions of 
standard materials or mathematical functions thereof, and corresponding properties of 
said materials, or products of said process ot' which .\ is a feed or yield of said process, 
the method comprises steps of 

25 measuring absorption at more tlian one wnvelength in the region 600-2600nm 

at the process line and producing absorption signals or mathematical functions thereof, 
by the spectrometer in accordance therewith; 

accessing the databanks of the computer in accordance with the absorption 
signals or functions thereof; 

30 comparing, by the computer, the absorption signals or functions thereof to the 

signals or functions thereof of the standard materials stored in the databanks, 

choosing at least one standard based on the comparing, said standard having 
the smallest average value of the absolute difference at each wavelength i between the 
signal (or function thereof) for the material and the signal (or flmction thereoO for the 

35 standards, with averaging of said properties or yields when more than one standard is 
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chosen; and 

controlling said process in accordance with the oiitputted property/yield. 

34. A computer programmed to perform the method ot'any of claims 1-32. 

35. Apparatus suitable for use in the method of any one of claims 1-33 which 
comprises an NIR spectrometer receiving at least one signal from a feed or product 
line in said process and being coupled to a computer to effect continuous measurement 
of the spectra of the feed and/or product and provide feed back control of the process. 
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